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STOCHASTIC HIGHER SPIN VERTEX MODELS ON THE LINE 


IVAN CORWIN AND LEONID PETROV 


Abstract. We introduce a four-parameter family of interacting particle systems on the line which can 
be diagonalized explicitly via a complete set of Bethe ansatz eigenfunctions, and which enjoy certain 
Markov dualities. Using this, for the systems started in step initial data we write down nested contour 
integral formulas for moments and Fredholm determinant formulas for Laplace-type transforms. Taking 
various choices or limits of parameters, this family degenerates to many of the known exactly solvable 
models in the Kardar-Parisi-Zhang universality class, as well as leads to many new examples of such 
models. In particular, ASEP, the stochastic six-vertex model, q-TASEP and various directed polymer 
models all arise in this manner. Our systems are constructed from stochastic versions of the R-matrix 
related to the six-vertex model. One of the key tools used here is the fusion of R-matrices and we provide 
a probabilistic proof of this procedure. 
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1. Introduction 

Integrable probability is an active area of research at the interface of probability/mathematical 
physics/statistical mechanics on the one hand, and representation theory/integrable systems on the 
other. Integrable probabilistic systems are broadly characterized by two properties: 

(1) It is possible to write down concise and exact formulas for expectations of a variety of interesting 
observables of the system. 

(2) Asymptotics of the system, observable and associated formulas provide access to exact descrip¬ 
tions of new phenomena as well as large universality classes. 

In light of these properties, there are two main goals in this area: 

(1) Build bridges between algebraic structures and probabilistic systems and in so doing, discover 
new integrable probabilistic systems and new tools by which to analyze them. 

(2) Study scaling limits of these integrable probabilistic systems and in so doing, expand and re¬ 
fine the scope of their associated universality classes and discover new asymptotic phenomena 
displayed by these systems. 
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In this paper we work to advance the first goal. We develop a four-parameter family of stochastic 
interacting particle systems which are built off of higher spin representations of the six-vertex model 
R-matrix. These systems benefit from properties inherited from the R-matrix. In particular, they are 
diagonalizable explicitly in terms of a complete set of Bethe ansatz eigenfunctions, and they also enjoy 
certain Markov dualities. We use these two facts to compute moment and then Laplace-type transform 
formulas for these processes. Asymptotics of these systems and the associated exact formulas are in line 
with the second goal defined above. We do not pursue this here, but note that such asymptotics have 
previously been performed on various degenerations of this family of systems [BC141IBCF1211BCR121 
IBCFV141ICSS141IOQ141 IBarl4l I Vet 141 IBCG14| . In all of those cases, the resulting phenomena were 
that of the Kardar-Parisi-Zhang (KPZ) universality class. It would be quite interesting to see whether 
other phenomena can be accessed beyond that of the KPZ class. 

Vertex models and more generally quantum integrable systems have long been objects of intense 
research within mathematics and physics (see, for example, the reviews |Fad96l lRes08j ). Generally, 
R-matrices are not stochastic, and nor are their associated transfer matrices. Here we work with an 
R-matrix which is stochastic (and arises from conjugating the associated transfer matrix). This enables 
us to define our vertex models on the entire line such as in |BCG14j . The associated transfer matrix is 
Bethe ansatz diagaonalizable, as follows either from taking a limit of the finite lattice algebraic Bethe 
ansatz or from the recent work IBorldl . Moreover, as opposed to on a finite lattice, there are relatively 
simple direct and inverse transforms with respect to these Bethe ansatz eigenfunctions 1BCPS14] (see 
also Appendix [Aj. Our proofs of duality are largely based on earlier methods from |BC13l ICorl4| . In 
both cases (diagonalization and duality) we first prove our results for the horizontal spin 1/2 (J = 1) 
and arbitrary vertical spin R-matrix. Then we construct higher horizontal spin R-matrices (J £ Z>i) 
via Markov functions theory, and the diagonalization and duality easily extends. This construction 
provides a probabilistic proof of the fusion of R-matrices on the line. Though initially the form of the 
R-matrix we find is not very explicit, we are able to demonstrate a recursion relation it satisfies, and 
then explicitly solve that in terms of g-Racah polynomials |Manl4j . 

We also have a more concrete motivation behind this work, which we now explain. The totally asym¬ 
metric simple exclusion process (TASEP) is arguably the paradigm for integrable stochastic interacting 
particle systems. The asymmetric simple exclusion procesqj (ASEP) and g-deformed totally asymmetric 
simple exclusion procesqj (g-TASEP) are both one-parameter generalizations of TASEFq In [BCS12] 
it was recognized that both systems enjoy Markov dualities and have moment formulas which can be 
written in terms of nested contour integrals. This prompted the question of whether both systems and 
their associated results can be united as special cases of a more general integrable probabilistic system. 
In jBCPS14] this question was addressed at a spectral level - it was shown that the eigenfunctions 
which diagonalize both systems are special cases of a more general set of eigenfunctions (recalled in Ap¬ 
pendix [A}. In this present work we provide a complete answer to this question in the affirmative. The 
four-parameter family of systems we introduce here have degeneration^ to both ASEP and g-TASEP, 
at the level of their transition operators. In other words, not only do the eigenfunctions degenerate, but 
so do the eigenvalues. 

In recent work, (CGRS14] discovered a family of higher spin particle systems which interpolate from 
ASEP (spin 1/2) to g-TASEP (spin infinity) and proved dualities for these systems using U q ( s^) sym¬ 
metry. This two-parameter family of system^ does not appear to be diagonalizable via Bethe ansatz 

^Particles jump left with rate (i.e., according to independent exponentially distributed waiting times of rate) p and 
right at rate q, assuming the destination is not already occupied. 

^Particles jump one to the right with rate given by 1 — g gap , where gap represents the number of holes before the next 
particle. 

^TASEP arises from ASEP by setting either p or q to zero, and from q -TASEP by setting q to zero. 

4 Saying that one system degenerates to another will mean that the second is accessed from the first either through a 
special choice of parameters, or through a limit transition. 

'’Spin and q being the two parameters. 
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KPZ equation / SHE / continuum polymer 


KPZ fixed point (e.g. Tracy-Widom distributions) 


Figure 1. Hierarchy of various degenerations of the four-parameter family of higher spin par¬ 
ticle systems. Each arrow represents a particular degeneration of models (with the exception 
of the arrow to the KPZ fixed point which is at the level of one-point limiting distributions). 

The arrows should be transitive, though not every additional arrow has been proved. 

except for the cases of ASEP and g-TASEP. Hence it is unclear whether many of the nice properties of 
ASEP and g-TASEP hold for the systems in between them. Besides the cases of ASEP and g-TASEP, 
the systems and dualities proved in [CGRS14] appear to be different than any degenerations of our 
four-parameter family of systems. 

Besides ASEP and g-TASEP, our family of particle systems has a number of other interesting degen¬ 
erations. Some of this hierarchy of degenerations is illustrated in Figure [T] wherein we show the relations 
to various known systems in the literature. There are, of course, many other degenerations to study. 
The exact nature of the degenerations are indicated below or in the cited literature. The left-hand 
side of the figure involves g-Hahn TASEP (PovlSl ICorl4j (see also Section 15.61) , discrete time g-TASEP 
[BC13j . g-TASEP [BG141lBGS12| . the strict/weak polymer |CSS141lQQ14| . and the semi-discrete Brow¬ 
nian polymer jOYOll lQ’C12j . The right-hand side of the figure involves the stochastic six-vertex model 
|BGG14| (see Section I5.5j) , ASEP [TW081ITW09IIBGS12] , and Brownian motions with skew reflection 
|SS14j . Both sides have limits to the KPZ equation / stochastic heat equation / continuum polymer 
IBG1 (ACQlll |MFRQ15| and yet further to the KPZ fixed point (e.g. with cube-root scaling and limit¬ 
ing GUE Tracy-Widom one-point fluctuations) |TW09[[AC^ITlIskTnliBGUl IBGF12UBGR.1 21 lBGFVT4l 
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ICSS141 100141 lBar!4l IVetl4l IBCGI4] , Besides these examples, there are also ‘determinantal’ particle 
systems (such as TASEP) which fall into the hierarchy of degenerations and are not depicted. 

Let us note that the log-gamma polymer |Sepl2[ ICOSZ141 fTLD14j . g-pushASEP |BP131 ICP15| . and 
two-sided g-Hahn ASEP [BC15| are not included in Figure [1] These systems are all diagonalized in the 
same basis (or degenerations thereof) as our four-parameter family of systems, and they have explicit 
and elementary eigenfunctions. While we expect that these systems can also be incorporated in some 
form into our hierarchy, we do not pursue this direction here and leave it for future work. 


1.1. Outline. Section [2] is devoted to the study of the J = 1 (or horizontal spin 1/2) vertex models. 
In particular, Section 12.11 provides definitions of the J = 1 i?-matrix, describes certain conditions on 
parameters under which the matrix is stochastic, and constructs three-parameter discrete time zero 
range and exclusion processes from the 17,-matrix. Section [2.21 provides the Bethe ansatz diagonalization 
of these processes. Section [2731 contains the statements and proofs of dualities enjoyed by these systems. 

Section [3] implements the fusion procedure through which we go from J = 1 to arbitrary J £ Z>i - 
thus yielding the fourth parameteio In particular, Section [3.11 explains the fusion procedure by which 
the horizontal spin is taken from 1/2 to J/2 for arbitrary J € Z>i. The higher spin 17-matrix is 
constructed through use of the theory of Markov functions which draws on certain special properties 
we check for the J = 1 case. The diagonalizability and dualities for the higher spin systems follow 
immediately from the J = 1 cases. The form of the 17-matrix (and hence higher spin zero range and 
exclusion processes) is initially rather inexplicit. Section [3.21 deduces a recursion relation in J satisfied 
by the 17-matrix. Section 13.31 provides an explicit solution to that recursion relation in terms of g-Racah 
polynomials (or terminating basic hypergeometric functions). 

Section |4] utilizes the duality between the J higher spin zero range and exclusion processes to com¬ 
pute nested contour integral moment formulas for the exclusion process with step initial data. These 
lead to a Fredholm determinant formula for an e q -Laplace transform which characterizes the exclusion 
process’s one-point marginal distribution. We do not pursue asymptotics, though this type of Fredholm 
determinant has been used before for such purposes (see references earlier in the introduction). 

Proposition 12.31 identifies four different cases of parameters under which the J = 1 17-nratrix is 
stochastic. Case (1) of that remark is assumed in the earlier sections of this paper. However, in 
Section [5] we explain how to extend our results to the other three cases. Section l5Tl addresses case (2). 
Section EJ identifies reflection and inversion symmetries of the higher spin 17-matrix. Section 15.31 
addresses case (3) and Section 15.41 addresses case (4). Section 15.51 demonstrates how the stochastic 
six-vertex model (BCG14) arises from our systems, and establishes a new self-duality of that model. 
Section 15.61 describes another degeneration to the g-Hahn processes |Povl3l ICorl4| . Section 15.71 briefly 
notes which results extend to spatially or temporally inhomogeneous parameters. 

Appendix [A] recalls key facts (namely, the Plancherel theory) about the Bethe ansatz eigenfunctions 
which diagonalize the systems considered herein. Appendix [B] contains explicit formulas for 17-matrix 
elements for J = 1,2,3. Appendix [C] states the Yang-Baxter equation satisfied by the 17-matrix. 

1.2. Acknowledgements. We thank Alexei Borodin for providing us with an early draft of [Bor El¬ 
and for early discussions related to this work as well as comments regarding an early version of the fusion 
proof described in Section 13.11 IC also appreciates discussions with Nicolai Reshetikhin regarding the 
general challenge of finding stochastic quantum integrable systems. IC was partially supported by the 
NSF through DMS-1208998, the Clay Mathematics Institute through a Clay Research Fellowship, the 
Institute Henri Poincare through the Poincare Chair, and the Packard Foundation through a Packard 
Fellowship for Science and Engineering. LP was partially supported by the University of Virginia 
through the EDF Fellowship. 


®The four parameters we work with are q, a, I, J or in a different parametrization, q, a, v, ft. 
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Figure 2. Left: The weight associated to a vertex with (clockwise, starting on the bottom) 

labels ii, Ji,*2, J2 in Z>o is given by Ra\i±, ji; Z2, j’2). The R-matrix takes as input (i\,j\) and 
produces output ( 12 ^ 2 )- Right: The representation of a vertex associated with (f 1 , J 1 : f 2 , . 72 ) = 
(3,1,3,1) in terms of arrows and particles. 


2. J = 1 HIGHER SPIN STOCHASTIC SIX-VERTEX MODEL 


2.1. Definitions and construction of processes. We will consider q, v as fixed throughout and thus 
only include other variables explicitly in our notation. We will use the notation Z>j = {n € Z : n > i}, 
Z<j = {n € Z : n < ?}, and 1# is the indicator function of an event Li. The symbol E will denote 
expectation with respect to the process or random variable that follows. 

We proceed now to define the R-matrix which will play a central role in all the follows. 


Definition 2.1 (The J = 1 R-matrix). For three generic complex parameters q, u, a and any four-tuple 
(*i, ji,*2 ; j2) € (Z>o) define a corresponding vertex weight as follows: For any m > 0, 


Ra\m, 0 ; m, 0 ) 

R^\m, l;m + 1,0) 


l + qg m 
1 + a 
1 - 

1 + a ’ 


R^\m,0;m — 1,1) 
R^{m, 1; m, 1) 


a(l — g m ) 

1 T o 

a + 

1 + a ’ 


( 2 . 1 ) 


and (ii, ji; * 2 > J 2 ) = 0 for all other values of (ii, ji, * 2 , 22 ) € (Z>o) 4 . Notice that all non-zero weights 
correspond to four-tuples such that ii+ji = i'l+j'i- a property we consider as being particle conservation. 
Weights are associated graphically with crosses labeled by (ii,ii, * 2 , J 2 ) in the manner of Figure [2] (see 
also Appendix [B]) . It will be convenient, at different points, to think about the i\-j \, * 2 , J 2 as recording 
the number of arrows (either up or right pointing) along the edges incident to a given vertex, or as 
recording the number of particles in a given location (the V s) and the number of particles to cross a 
given edge (the f s). 

We will treat these weights as matrix elements. Define a vector space V 1 with basis elements 
{0,1,...,/} if I € Z>i, and {0,1,...} otherwise. Likewise define a vector space ]H J with basis ele¬ 
ments {0,1,...,/} if J € Z>i, and {0,1,...} otherwise. As a matter of convention, when describing 
a linear operator acting between these spaces (or their tensor products) we will only describe matrix 
elements in the above basis. Then for I e C such that u = and J = 1, Ra^ : V 1 <g> H" 7 -A Y 1 <S> 
is defined by its matrix elements R^\ii,ji; / 1 , J 2 )- We will associate /?-matrices to vertices (x,y) € Z 2 
and denote vector spaces [Y J ] X and [H J ] 2/ as associated with column x and row y. We write [Ra^]x,y 

to mean the matrix which acts as Ra^ on \V J ] X <8> and the identity on all other [V 7 ]^ and [H 17 ]^. 

This convention is illustrated in Figure [3l 
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[M J ] 

[H J ] 

[H J ] 

[U J ] 



Figure 3. The association of //-matrices to each vertex in Z 2 as in Definition O 

Remark 2.2. These weights are closely related to the matrix entries of the horizontal spir@ 1/2 (i.e. 
J = 1) six-vertex model //-matrix. For instance, changing variables to s,u via a = —su and v = s 2 , we 
can match the above defined weights with those of [Borl41 Definition 2.1]: 

Ra\h,ji;i 2 ,j 2 ) = w u (i 2 ,j 2 ;ii,ji)(-s) : ' 1 (-suy 2 ~ 31 . 

The factor (— su) 32 ~ n is a conjugation which does not affect the overall transfer matrix (cf. Remark l2.8D 
whereas the factor {—s) n corresponds to conjugation of the transfer matrix (cf. Remark 12.71) . [Borl4l 
Proposition 2.4] further matches w u (i\, ji; J 2 ) to the matrix entries of a particular normalization of 
the //-matrix considered in |Manl4] when J = 1. 

We now identify various ranges of parameters under which our ill-matrix is stochastic. Notice that 
our weights have been normalized so that for ii,j\ fixed, 

111*2,32) = 1 - 

*2U2 

If q. v, a are chosen so that all weights are non-negative, then Ra' 1 is a stochastic matrix and provides 
the transition probabilities for going from the pair (ii,ji) to (* 2 , J 2 ) (see Figure [2] for a graphical 
representation of this transition from inputs to outputs). 

Proposition 2.3. The following choices of parameters ensure non-negativity (and hence its stochastic- 
ity) of : 

(1) q, v E [0, 1), and a > 0, 

(2) q € (-1,0], a G (0,l/|g|), and v € ( - l/\q\, min(l, a/\q\)), 

(3) q G [0, 1), v = qfor I G Z>i, and a < —q" 1 , 

(4) q G (1, + 00 ), v = q~ I for I G Z>i, and —q~ I < a < 0. 

Proof. Each case follows from straightforward inspection. □ 

Tn general, for I, J € Z>i one says that our i?-matrix is vertical spin 7/2 and horizontal spin J/ 2. 
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Left to Right Update 
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Figure 4. A schematic showing the discrete time zero range process sequential update. The 
grey particles on the bottom of the figure are the ones which will move to the right by one. 


Remark 2.4. In the third and fourth cases, the weights are not always non-negative, however restricted 
to i G {0,...,/} (i.e. Y 1 ) they are. The only way to transition out of this range for i is to utilize the 
weight Ra\l,l;I + 1,0), but due to our choice of v, this is zero. The first choice relates to models 
generalizing the g-Boson stochastic particle system |SW98l IBCS12] whereas the third and the fourth 
examples relate to models generalizing ASEP and the stochastic six-vertex model (coming from letting 
1 = 1) )BCG14| . 

We will assume that q } v,a satisfy case (1) above, namely that 

9 , 1/6 [ 0 , 1 ), and a > 0 , ( 2 - 2 ) 

and prove all of our results under those conditions on parameters. In Section [5] we describe how our 
various results extend to the other choices of parameters from Proposition 12.31 

We turn now to define certain sequential update, discrete time zero range processes and exclusion 
processes based off of our stochastic R-matrho This is a general construction whose only input is the 
stochastic R-matrix. We will make use of this construction later when Rp^ is replaced by R^'P . 

Definition 2.5 (State spaces). Define the space of left-finite particle configurations on the line as 

G = j<? = (■ ■ ■ , g-i, go, gi, ■ ■ •) : all gi € Z>oU{+oo}, and there exists iGZ such that g t = 0 for all i < x 

and likewise the space of right-finite particle configurations on the line as 

Y = |y = (• • • , y~i, yo, yi, • • ■) : all yi € Z>oU{+oo}, and there exists x 6 Z such that yt = 0 for all i > x 

^Strictly speaking, these might not be the most accurate names for these processes, but for lack of a better name they 
will suffice. 
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For k E Z>i, define the space of £>particle configurations on the line as 

= | g= (••• ,g-i,g 0 ,gi ,...) : ^.g* = A:}, and Y k = { y = (• • • , y-i, yo, yi, ■ ■ ■) : = &}■ 

i i 

Though G k and Y k are the same, we differentiate to keep track of the separate Markov processes for 
which they will be the state spaces. The spaces G k and Y k are (respectively) in bijections with (W for 
Weyl chamber) 

WG fc = | rh = (mi < ■ ■ ■ < mrik) : all nii E z|, and WY fc = |ft = {n\ > ■ ■ ■ > rik) '■ all rtj € z|. 

The bijection is given by associating to a state g or y, the ordered list of the k particle locations. As a 
convention, for g we associate rh with weakly increasing particle location order, and for y we associated 
n with weakly decreasing particle location order. Given an operator B acting on functions from G k to 
C, we will overload notation and let B also denote the operator on functions / from WG k to C defined 
via ( Bf)(m ) = ( Bf)(g ) where g and rh are associated via the bijection, and f(g) = f(rh). Finally, 
define (for later use) Gj,Yj,G k , and Y k to be restrictions of the respective spaces so that each g j or yi 
lies in {0, ...,/}. Define WG k and WY k as the respective images of G k and Y k . In other words, having 
no clusters of equal mi’s or m 1 s of length large than I. 

Define the space of right-finite exclusion particle configurations 

X = jx = (xi > X 2 > ■ ■ ■) : all x, E zj, 
and for > 1 define the space of A r -particle configurations 

X N = jx = (xi > • • • > xn) : all Xj E zj. 

By convention, we define ‘virtual particles’ Xi = +oo for all i € Z<o- 

Definition 2.6 (J = 1 high spin zero range process). For k > 1, define the /c-particle discrete time 
J = 1 higher spin zero range process g{t) with state space G k according to the following update rule 
(see Figured]). Given state g we update to state <f sequentially.. Start at the left-most site x € Z such 
that g x > 0. Let h x = 0 and randomly choose g' x and h x+ \ according to the probability distribution 
Ra\g x ,h x ] g' x ,h x+ i). Now proceed sequentially so that given g x +i and the choice of h x+ \ randomly 
choose g' x+ i and h x+ 2 according to the probability distribution Ra\g x + 1 , h x+ i] g x+1 , h x+ 2 ). Continue in 
this manner, augmenting x. Since there are finitely many particles and since all probabilities of the form 
-Rq)(0, ji; iii J 2 ) are strictly less than 1, it follows that eventually the output of the sequential update 
step will be all zeros. The update can be stopped and all subsequent g' values are set to zero. We 
will denote the transition probability (given by the product of /?,-matrix weights) from state g to (f as 
J3 a,qa (g, if) and the transition operatoijj with matrix entries as B a,qa . This operator has a well-defined 
action on functions / which are bounded as g goes to infinity!^. The odd superscript notation B a,qa 
is due to our eventual extension to B a ^ for arbitrary j3 (cf. Remark 13.181) . Note that the dynamics 
of this process preserves the total number of particles due to the particle conservation property of the 
/I’-inatrix. 

We also define the infinite-particle version of the process g(t) with state space G. Since the update 
is from left to right and we are dealing with left-finite initial data, we can show that the update is well- 
defined. Towards this end, for M E Z, define the restriction of a state g to Z<m as g\ M = ( gi ■ lj<iw)iez 
and define a sequence of finite-particle processes g\ M (t ) such that g | M (0) = g\ M - Each of these initial 
data has finitely many particles, and hence the evolution according to the above defined finite-particle 
version of the J = 1 higher spin zero range process is well-defined. We define g{t) as the inverse limits 
(in law) of the restriction of g\ M (t) to Z<m (Note: even though the initial data was restricted to Z<m, 

9 The operator acts on functions / : G k —> C as (B a ’ qa f)(g) = B a ' qa (g, g')f{g') 

^ 9 By g going to infinity, we mean that the right-most particle in g goes to infinity. 





4 i ) ( 3 , 0 ; 2 , 1 ) ^>( 1 , 1 ; 2 , 0 ) ^‘>( 2 , 1 ; 2 , 1 ) 

--sA ry^ Wry^ ---sA ry^ 

/ T«Ay 3 ^ ^^2 


^4 X3 X2 


Right to Left Update 

Figure 5. The exclusion process constructed from the R-matrix. The first particle x\ is 
updated and then, based on the length of the previous update, each subsequent particle is 
updated. 

the evolution may have left that sector). Because of the left to right update and one-sided nature 
to particle movement, these restrictions form a consistent family and the desired inverse limit exists. 
Though we cannot write transition probabilities, we can still define a transition operator acting on a 
suitable domain of functions. The operator B a,qa will act on functions / : G —» C which are stable 
at infinity. By stable at infinity we mean that for each g E G, over all g 1 > g, f((f) is uniformly 
bounded and /(<f| M ) converges uniformly to /({?)■ The inequality g 1 > g means that for all x E Z, 
Yli<x 9i — Yli<x 9i (i n other words, the state g 1 can be accessed from g via moving particles to the 
right). For such stable functions, liniA/^-oo (B a,qa f) {g\ M ) exists for all g, and defines ( B a,qa f)(g ). 

In the same manner as above, we define space reversed zero range processes y{t) in which particles 
are updated right to left with state space Y k (or ¥). This reversed process involves the R-matrix 
which is reflected in the y-axis. The space reversed transition operator B a,qa = PB a ’ qa P~ 1 where 
(R/) ((yi)iez) = /((y- i)i G z) is the space reversal operator (note that P 1 = P). We may likewise 
extend from finite to infinite particle configurations. 

Remark 2.7. The J = 1 higher spin zero range process describe in Definition 12.61 should be thought 
of as a full line version of the transfer matrix built from the R-matrix with matrix elements given in 
Definition 12.11 The standard construction of a transfer matrix on Z/LZ involves taking the product 
of R-matrices and then tracing out the horizontal space. In other words, one defines the transfer 
matrix as trpj^ ([Rq ^]i,i[Ro ■ ■ • [Ra . The resulting matrix maps [V^i <8> <8> ■ ■ ■ <8> [V 7 ]q to 

itself. Under standard normalization, it is not clear how to directly construct transfer matrices on Z. 
However, if we work with stochastic R-matrices, then all weights are strictly less than 1, except for the 
weight Rq 1 ) (0,0; 0,0) = 1. This enables us to make sense of the infinite product of these (stochastic) 
R-matrices, at least when restricted to the sector in which there is a finite total number of particles. As 
in Definition 12.61 with a little more work one can likewise construct infinite particle number versions of 
these resulting stochastic transfer matrices (see also |BCG14] L 

Remark 2.8. The construction of the finite particle zero range process given in Definition 12.61 is 
invariant under conjugation of the R-matrix elements via multiplication by /(ji) j for any non¬ 
zero function / : BU —>• C. Such a conjugation may, however, destroy the stochasticity of the individual 
R-matrices (despite retaining that of the entire transfer matrix). 

Definition 2.9 (J = 1 higher spin exclusion process). Via a gap/particle transforirF*!. we define an 
exclusion process. For N > 1, define the V-particle discrete time J = 1 higher spin exclusion process 
x(t) with state space ~K N according to the following update rule (see Figure [5]) . Given state x we 
update to state x! sequentially. Start at x\, randomly choose x\ E {xi,xi + 1,...} according to the 

may consider our zero range process as describing the size of gaps between labeled particles. This equivalence is 
unique up to an overall shift in the labeled particle positions. 


R^(oo,0; oo, 1) 



X\ 
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probability distributioi£3 Ra\oo, 0; oo, x[ — aq). Proceed sequentially up to j = N so that given 
x'j and Xj + \. randomly choose x' +1 E {xj + \,... ,Xj — 1} according to the probability distribution 

Ra\xj — Xj -|_i — 1, x’j — Xj ; x' 3 — x'j +l — 1, x' j+1 — x j+ i). We will denote the transition probability from 
state x to x' as T a ’ qa (x,x'), and the associated transition operator with these matrix entries as T a,9 “. 
This transition operator has a well-defined action on functions which are bounded as x goes to infinity. 
In order to define the infinite version of this process with state space X it suffices to observe that the 
first N particles evolve according to the IV-particle exclusion process (and hence the infinite-particle 
version can be defined as an inverse limit of consistent laws of these processes). We can also define the 
transition operator’s action on functions / : X —>• C which are stable at infinity. By stable at infinity we 
mean that for each x E X, over all x 1 > x , f(x') is uniformly bounded and f{x'\ N ) converges uniformly 
to f(x'). The inequality x! > x means that for all i E Z>i, x\ > Xi, and x'\ N = (aq, ..., aqv)- For such 
stable functions lim 7 v_ 5 . 00 (T a,qa f) (x\ N ) exists for all x. and defines (T a,qa f)(x). 

Definition 2.10 (g-Hahn distribution and zero range process). For generic complex q , y, u, and y E Z>o, 
define the (complex) probability distributiorFI on s E {0,... , y} as 


y) = Vq,^A s \y) = R 


{u/y 1 q)s(y 1 q)ys {q\q) j 


Aq)y 


(q-,q)s{q;q) y -s 


with {a-,q)n := JJ(1 - aq j x ). 
3 = 1 


For k > 1 fixed and i E Z, define linear operators [H^\ ■ via their action on functions / : Y fc —>• C 


Hi 

([u a ]j)(y) = v-A s i\yi)f{y^i~i)- 


Si= 0 


Note the inclusion of the negative sign in (f- a {si\yi)- Here, for y E Y k , we have set yf = (..., + 

s,yi — s,... Define the fc-particle space reversed q-Hahn zero range process transition operatoiFl 

n a = ..- [H a ]_ 1 [H a ] 0 [H a \ 1 --- . 

In the same way as for B a,qa (up to space reversal), we also define the transition operator’s action 
on functions / : Y —>• C which are stable at negative infinity. Likewise, define the non-space reversed 
operator "H a = P'H a P ~ 1 where P is the space reversal operator from Definition 12.61 

Remark 2.11. We will use the following properties of fp^(s\y) which one can readily check (here 

s € {!,... ,y}): 


t i \ , f i \ i- q y+1 s l~vq s 1 /n 0 ^ 

<pA a \y) = w - !\y) ■» l _ qS — i-nqv-s » ( 2 - 3 ) 

vqA s \y) = (C 1 - w v - a )<pA s \y) ~ A 1 -q y+l ~ s )vA s -My))- ( 2 - 4 ) 


12 This corresponds to taking the well-defined limit m —>■ oo in Definition ^. II 

1: ^One condition under which this is a bona-fide positive probability distribution is if |g| < 1 and 0 < v < y < 1. For 
discussion and references regarding this distribution, see |Povl3l ICorl4| . 

*^If the parameters q, y, v are such that is always positive, then [H 01 ] corresponds with the Markov update by which 
Si particles are moved from position i to i — 1 according to the distribution <p- a (si\yi), and T-L a is the Markov update 
by which one updates each site in parallel according to the H a single-site update. We do not, however, rely upon such 
positivity in our use of this distribution and its associated operators. 
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2.2. Bethe ansatz diagonalization. In Appendix [A] we recall the Bethe ansatz eigenfunction which, 
according to the result we now prove, diagonalize the higher spin zero range process with transition 
matrix B a,qa . We will rely upon the Plancherel theory developed in [BCPS14] for these eigenfunctions 
(also reviewed in the appendix). In what follows we will, as described in Definition ^. 51 overload notation 
and let operators acting on functions of the y variables also act on functions of the n variables via their 
bijective association. 

Proposition 2.12. Assuming | jEff. f+^| < 1 f or 1 < * < k, theift^ 


{ ' g)W II 1 + 

i— 1 


Proof. We appeal to the known eigenfunction relations for the J = 1 higher spin six-vertex model. This 
can be derived on the periodic lattice via algebraic Bethe ansatz (cf. [ Res08l 1. We use jBor!4l Corollary 
4.5 (i)] wherein it is shown (via a symmetric function theory approach) that if 


Z{ — s u — s 
1 — SZi 1 — su 

for alii E {1,..., k} , then 


n 


1 — quzi 
1 — UZ{ 


c (n)Fft(zi ,..., z k ) = ^2 G h , /H ( v)c(n')F H ,(zi,..., z k ). 

n' 


Here ft 1 , n E WY k and, after changing z K > sz , 


c(n)F H (szi,...,sz k ) 


id 

i= 1 


— S 


\rii 


- S 2 Zi 




where v = s 2 in Th(n). According to |Borl4l Definition 3.2] the term G^r/^(v) is equal to the product 
of the 17-matrix weights from n' to n where the weights are given as in |Borl4l Definition 2.1]. Remark 
12.21 explains the relation of the weights in |Borl4j to those considered herein and setting a = —sv and 
u = s 2 yields 


i= 1 


+ azi 


J g Q '9 “) r \If 


n 


This, likewise, implies the desired relationship for left eigenfunctions as well. 


□ 


Proposition 2.13. For a E C and 1 — vzj / 0 for 1 < j < k, 

(«“*!)(«)=nP-^4(«) 

i =i 1 Zl 

Proof. This follows from |Povl3i fCoiT4l IBCPS14] , In particular, [BCPS14] Proposition 5.13] records 
the desired result for 0 < u < —a < 1. The operator PL 01 depends polynomially on ct, as does the 
eigenvalue. Thus, since both sides above are polynomial in a and equal for an interval of values, they 
must match for all a E C. □ 

In Section 15.61 we develop the relationship between the operator PL a and higher spin versions of the 
zero range process transition operator. 

In order to understand the following corollary, the reader is encouraged to recall from Appendix lAl 
the direct and inverse transforms F q ’ v , J q,v and the space W^ax on which J q ' v T q,v acts as the identity. 


l J By applying the space reversal operator P we can also produce eigenrelations for B a ’ qa . 
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Corollary 2.14. On the space W k , aY , B a,qa = (%“) l PL qa , B a,qa = ("%“) 1 'H qa , and 7i a commutes 
with itself for different values of a. 

Proof. This follows by spectral considerations. Let’s prove the first claim. The formula for J q,u implies 

that@ 

((«")- 1 /)(s) = /)(■)) (a) 

(«’”/)(«) = ■7"'(e^(0(J™'‘7)(-))(»), 

where ev a (z) = n!=i i-^z* • Combining these yields 

( (y. a )~ 1 'H ga f^j (n) = J q ’ u (ev a , qa (-){^f)(-))(n) 
where ev a ,b(z) = IIi=i l+a^z- • On the °tber hand, it follows from Proposition 12.121 that 

(B a ’ qa f)(n) = J^(ev a , qa (-){^f)(-))(n) 

as well. Conjugating everything by the space reversal operator P produces the second claimed result. 
Similar considerations and Proposition 12.131 imply the last commutation relation. □ 

Remark 2.15. One of the methods used in Bethe ansatz is to rewrite an operator as the direct sum of 
one-dimensional operators subject to two-body boundary conditions. Not every higher spin transition 
operator is amenable to this method - for instance, the stochastic six-vertex model }BCG14j . This 
proposition shows that B a,qa (and eventually via fusion, the general B a ’@ operator) can be written as 
the ratio of g-Hahn transition operators. It would be interesting to see if this fact goes through to the 
case of transfer matrices on the finite lattice Z/LZ. 

2.3. Self duality. 

Definition 2.16. We define a number of duality functionals. For x € X and y E Y k for some k E Z>i, 
defirnO 

iez 

For g E G and y E Y define 

G(g, y) = q^i>j 9iVj , and G (g,y) = 9iy P 

Notice that unless g E G k and y E Y k ' for some k, k! E Z>i, G (g,y) will equal zero. In the case that 
g € G k and y E Y k ' notice that 

G(<7,y) = G(g,y) ■ q kk '. (2.5) 

For g E G k and y E Y k ' we may overload these functionals as described in Definition 12.51 by replacing g 
by m E WG fc and y by n E WY fc . For m£Z (or n E Z) and g E G (or y E Y) 

N m(<?) = ^2 9e, and l\lj n (<?) = ^ ge- 

l<m t>m 

For g E G k bijectively equivalent to m E WG fc and y E Y we have 

k ^ k 

G(y,y) = G(m,y) = ]Jg Nm i- l(?/) , and G(g, y) = G(m, y) = q~ N ™i( y) . 

i =1 i =1 

the below formulas, the dot represents the variable integrated in the application of the inverse transform. 
l^We employ the convention that the product is zero if yi > 0 for any i < 0 . This is in accordance with the convention 
that Xi(-) = +oo for i < 0 . 
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For y E G and y E Y fc/ bijectively equivalent to n E WY fc we have 


G (y, y) = G(y, n) = g N "i+ l(9) , and G(y, y) = G(y, n) = y 


-<(£) 


2=1 


2=1 


Recall the (/-integer [n] = . For fh E WG 1 ’ and y E Y dehne 


G m (m,y) = ]^[[y mi ] 9 • y N "u- l(2/) , and G m (rh,y) = JJ[; 


,2/mJq • Q 


-N L(?) 


2=1 


2=1 


and for y E G and ft E WY fc define 


fc' 


G n (</,«) = n^- 1 ], • <A i(9) , and G n (y,n) = JJ[y ni ] 9 • q 


-<(s) 


2=1 


2=1 


Definition 2.17. Call x E X well-adapted to H if for y E Y fc the function y t —> H(x, y) lies in the space 
W^ax (see Dehnition lA.5l) . Call y E G well-adapted to G if for y E Y fc the function y i->- G(y, y) lies in the 
space W^ ax . Call y E Y well-adapted to G if for g E G fc the space reversal of the function g i—)■ G (y, y) 
lies in the space 


Remark 2.18. There are a number of readily accessible examples of well-adapted x, g. and y. Step 
initial data z* = —i, i E Z>i is well-adapted to H. To show this, note that for c > 0 small enough 
and C > 0 large enough the function y H(x, y) lies in WY, p(c C) which, by Corollary IA.6I is a subset 

of W^ ax . Similarly, one shows that spiked initial data y^ = +ool,=o and likewise y* = +ooli=o are 
well-adapted to G. 


We now state and prove our first duality result, one between the exclusion process x(t) and zero range 
process y(t). The approach of the proof follows, for the most part, that of the proof of the discrete 
time Bernoulli y-TASEP duality in [BC13] . Indeed, setting u = 0 our exclusion process becomes the 
Bernoulli y-TASEP. 


Theorem 2.19. The J = 1 higher spin exclusion process with transition operator T a ^ a is dual to 
the space reversed J = 1 higher spin zero range process with transition operator B a ' qot with respect to 
H (x,y). Precisely, 

T a ’ qa H = H (B a ’ qa ) T , 

where the equality holds for all matrix elements indexed by x EX well-adapted to H and all y E Y k . 


Proof. We will demonstrate the following ‘implicit’ version of the self-duality. Consider the exclusion 
process x(-) from Definition 12.91 with some (possibly random) initial data x(0) E X and define 


It(y) = E H(x(t), y) 


where the expectation is over the evolution of x(-) (as well as the initial data, if it is random), and 
y E Y k for some k E Z>i. Then, assuming the initial data is such that R(y) is everywhere finite, we 
claim that following equality 


{n a i t+ i)m = ( n qa i t )(y), 


( 2 . 6 ) 


holds for all t E Z>o and all y E Y fc . 

Before proving (12.6p . let us observe how it readily leads to the proof of the theorem. In matrix 
notation, (12.61) implies that 

T“’ 9Q H(R“) t = H {U qa ) T , 

where we note that H is easily seen to be in the domain of these transition operators. Since in the 
hypotheses of the theorem we have assumed that x E X well-adapted to H, we can appeal to Corollary 
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12.141 (which shows that B a,qa = ("H“) 1 'H qa when acting on function in W^x) to complete the proof 
of the theorem. 

Thus, it remains to demonstrate (I2.6P and the remaining portion of this proof is devoted to that goal. 


Let N be the smallest integer such that all y% = 0 for i > N. Note that by the convention on H 
being zero if yi > 0 for any i < 0, it suffices to restrict the product i < N to 0 < i < N and modify the 
definition of T~L a to only include terms [%“]. for 1 < i < N. In other words, we can treat 0 as a sink 
since our functional becomes 0 for all i < 0. We may now argue similarly to |BC13| . 

For j > 1 and t > 0, let a J t denote the sigma-algebra generated by the random variables x\ (t), ..., Xj(t). 
Conditioning on the history of the whole process up to time t, we will show that, as cr^ measurable 
random variables, 

H a E 

Since operators TV* and T~L qa in both sides have a sequential structure (corresponding to first moving 
si particles from y\ to yo, then s 2 particles from 7/2 to y \, etc.), we can further condition on what 
happened to particles aq,..., during time step t —>• t + 1 (for any i = 1,... , N). We will show that 
the relevant contributions to both sides of (12.71) behave as they should (i.e., the parameter a in the 
operator at time t T 1 is replaced by qa at time t). That is, we will show that for i > 2, 


N 

n- 

. 1=0 


(xi(t+l)-M)j/i | N 


N 


n qa n q(xi(t)+i)yi 


(2.7) 


i =0 


Vi 

Si=0 


q(xi(t+l)+i)(yi—Si)q(xi—i(t+l)-\-i—l)si 


a 


N 
t > 


a 


i— 1 
t +1 


Vi 


— 22 V-qa 

Si =0 


I Vi )q ( ' Xi ^ Vi - s i) ; 


and for i = 1, 


( 2 . 8 ) 


p_ a (0|i/i)E 


7 (xi(t+l)+l)i/i I a N 


= V- qa (0\yJq^W+Vn. 


(2.9) 


First, note that (12.91) is straightforward: conditioned on the knowledge of x\ (t), the first particle jumps 
to the right by one with probability a/(l + a) and stays put with probability 1/(1 + a). Therefore, 


<P- a (0\yi)E 


(xi(t+l)+l)3/i I N 


= v?_ Q (0|yi) 


a 


1 T cr 


q(xi(t+i)+i)yiqVi _)_ 1 g(xi(t+i)+i)yi j 


1 + a 




= <P- qa (0\ yi )q^V +1 ^. 


Here we have used s = 0 case of (12.41) with y = —a, which reads 

V- qa (0\y) = —(1 + aq y - s )ip- Q (0\y). 

To show (12.81) . denote by I the indicator of the event that Xj_i(t + 1) = Xj_i(t) + 1, i.e., that the 
particle Xi-i has jumped to the right by one during time step t —» t + 1. This indicator is 
measurable, and it will help us to compute the conditional expectation in the left-hand side of (12.81) . 
Also for any r,s > 0 denote Z rs := q{ x i(i-)+ l ) r q( x i- i(t)+*-i)^ Using the definition of the dynamics of x, 
we can write 


E 


q(xi(t+l)+i)r q(xi-\(t+l)+i-l)s 


(J 


N 
t 1 


a 


2—1 

t +1 


= Iq S Z T)S 


r a + vq Xi ~ 1 1 _ jyqXi-iitf-xiit)-!' 


IT OL 


IT ol 
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f r a(l — 1 + ag *i-i(t)-*i(t)-i\ 

+ (1 — l)/r }S q -——-1-——- . 

\ 1 + a 1 + a / 

Noting that q x i-i( t )~ x i( t )- 1 Z r , s = Z r _i. (S+ i, we can simplify the right-hand side above to 
= -f- (W+ s + - aq r - l)Z r , a + ( vq r+s - V + aq r - a)Z r _ ljS+1 ) 

+ Y + a (( ac l' + 1 )Zr,s + i~ a Q r + a )Zr—l,s+lJ 

(note that when r = 0, the coefficient by Z r - i, s +i is zero.) The left-hand side of (|2.8|> is then equal 
to the sum (over Sj) of tp- a (si\yi) times the above expressions with r = yi — Si and s = s*. That is, 
left-hand side of (12.811 takes the form 


I 


Vi 


— y z, 

1 + a ‘ 


Ui-Si,Si 


Si =o 


Vi 


+ 


1 H - Q( 


E z < 


2/i Si,Si 


ip- a (si\yi)(aq Vi + q Si - aq Vi Si - 1) 

+ ip- a (si - 1| Vi){vq yi ~ vq Si ~ l + aq Vi ~ Si+l - a) 
V-a(si\yi)(aq yi ~ Si + 1) + <f- Q (si - 1| yi)(-aq Vi ~ Si+1 + a) 


Si=0 


By properties of <f- a (|2.3D . (12.41) . the expression in the square brackets in the first sum vanishes for 
any sf, and the expression in the square brackets in the second sum is equal to <P-qa(si\yi)■ This yields 


Having now established (12.81) and (12.91) . we can now prove (12.71) . We have (assuming y 0 = 0, otherwise 
(|2.71) is trivial) 


[«“]i...[W“] JV _i[W“]„E 


(xi(t+l) + l)j/l]g 


(x2(t+l) + l)3/2 


,E 




1 

N-1 

1 

to 

1 

1 

1 

1 a t 

a t+ 1 

1 a t 

°t+l 

... | <J t 

a t+ 1 

\ a t 


= [W«]i...[W«] s _ 1 E 


(xi(t+l) + l)l/ljg (X2(t+l) + l)y2 


N -1 (t+ 1)+ N— 1) VN -1 


VN 


y C P-a(s N \y N )E 

sn=0 

= [«“], •••[*“] 


^(*jv(t+l)+JV)(s/jv-sjv)g(®i\r-i(t+l)+JV-l)sAr | „JV-1 


1 

N-i 

1 

CM 

1 

1 

1 

1 

a t 

a t +1 

a t 

°t +1 

... | <T t 

a t +1 

a t 


(®i(t+l)+l)j/i]g 


(®2 (*+!) + 1)»/2 


...E 

Vn 


Q 


1 N 

to 

1 N 

l 

1 JV 

1 a t 

a t+1 

... \ a t 

a t+ 1 

1 °t 


X ^ ^_ gQ ( SJ v|y7v)g (a;jvW+JV)(2/iV “ Sjv) Q (3:jv - lW+JV_1)sjv . 

•sjv=0 

The first equality is by definition, and the second equality is by an application of (12.81) corresponding 
to i = N, which leads to replacement of the operator [?7 Q: ] (for time t + 1) by the operator [H qa ] N 
(for time t). Continuing using (12.81) for i = N — 1,..., 2 and (12.91) for i = 1, we arrive at the desired 
identity (12.71) . and hence complete the proof of the theorem. □ 

We turn now to self-dualities of the zero range process. The proofs are considerably less involved and 
rely on an earlier discovered identity [Corl4i Proposition 1.2]. 
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Theorem 2.20. The J = 1 higher spin zero range process with transition operator B a,qa is dual to the 
space reversed process with generator B a,qa with respect to G (g,y) as well as G (g,y). Precisely, 

B a ’ qa G = G (B a ’ qa ) T , 

where the equality holds for all matrix elements indexed by g E G k and y E Y, or by g E G and y E Y k . 
Likewise, 

B a ' qa G = G (B a ’ qa ) T , 

where the equality holds for all matrix elements indexed by g E G and Y, provided both sides of the 
above equation are finite. 


Proof. We prove this theorem in a few stages. Initially we deduce an ‘implicit’ form of the G duality 
in (I2.1UI) . From this we deduce the G duality claimed in the theorem. Finally, we use the G duality to 
deduce the G duality. 

We begin by proving the following implicit form of the G duality. For any k E Z>i, 

U qa G(U a ) T = H a G(H qa ) T , (2.10) 

where the equality holds for all matrix elements with g E G k and y E Y well-adapted to G, or with 
g E G well-adapted to G and y E Y k . We will assume the second case below, though the first case 
follows similarly. 

To show (12.101) . we demonstrate first that 

G{U^) T = U^G ( 2 . 11 ) 

for all g € G and y E Y. This follows from the fact that for each i E Z, 

Vi Vi 

G([H li ]i) T {g,y) = ip /1 {s i \y i )q Si9i G(g,y) 

Si= 0 Si=0 

9i 9i 

= ^2 Tu.(ri\gi)q riVi G{g,y ) = ^ ^(ri|ft)G(^t +1 , y) 

n=0 r;=0 

= m t G(g,y). 

In the first and last equation above we have written the composition of operators followed by (<?, y) to 
denote the corresponding matrix element. The equality between the end of the first and beginning of 
the second line above relies on an identity proved in |Corl4l Proposition 1.2]. 

From (12.111) . we find that 

'H qa G('U a ) T = G (H qa ) T {H a ) T = G {U a ) T (U qa ) T = H a G (U qa ) T . 

The only step in this deduction which requires justification is the commutation relation (T~L qa ) T ('H“) T = 
('H a ) T ('H qa ) T . This, however, follows from Corollary 12. 141 and the fact that we have assumed that g£G 
is well-adapted to G and y E WY k . Thus, we have established (12.101) . 

Under the assumption that g € G is well-adapted to G and y € Y k we may apply the first two 
identities from Corollary 12.141 to equation (12.101) to deduce the G duality statement in the theorem. 
The theorem asks for this duality to hold for all g € G (without the well-adapted condition). Indeed, 
all g € G with = 00 are well-adapted to G. Thus, it remains to show that we can extend the 

duality to g E G k for any k! E Z>i. For g £ G k fixed, and M sufficiently negative (so as to be less 
than the location of the most negative particle in g) let g AI = +ool m + 9, where 1 m is the vector of 
all zeros, except a one at M. In other words, g AI is equivalent to g except with an infinite number of 
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particles added at site M. From the above argument we know that the G duality in the theorem holds 
for matrix elements g M £ G and any y £ Y fc . It remains to show that both 

lim B a ’ qa G(g M ,y)=B a ’ qa G(g,y), 

M—>—oo 

lim G(B a,qa ) T (g M ,y) = G(B a ’ qa ) T (g,y). 

M—t—o o 

Let us justify the first limit, as the second follows similarly. Call f(g M ) = G(g AI ,y) and recall that 
(B a,qa fMg M ) gives the expectation of / after one step of the zero range process started from initial 
data g AI . The only way that the infinite number of particles at M can affect the value of / is if one of 
them makes its way past the left-most particle in y. However, this requires a large number of (0,1; 0,1)- 
vertices. Since the weight of these vertices is strictly less than one, this probability goes to zero as 
M —>• —oo. Since the above defined / is bounded by one, the desired convergence result clearly holds. 
This establishes the G duality in the theorem. 

Turning to the G duality, let us first prove it for g £ G k and y £ Y k ' for some k, k' € Z>i. Recalling 
(12.51) we have that G(g,y) = q~ kk G(g,y). Since q~ kk is a constant, multiplying the G duality by it 
yields the G duality. We now extend to all g £ G and y £ Y such that both sides of the G duality are 
finite. Under these conditions there must be a finite number (say k) of particles in g which lie to the 
left of some particle in y, and likewise a finite number (say k') of particles in y which lie to the right 
of some particle in g. It is easy to setpj that replacing g with its k left-most particles and replacing y 
with its k' right-most particles, the value of the left-hand and right-hand sides of the G duality identity 
are unchanged. This completes the proof of the duality and hence the theorem. □ 

Remark 2.21. The duality involving H and G can be made to look rather similar, though, to our 
understanding, they are not equivalent. From a state x £ X define its gaps via eg = Xj_i — Xi — 1. Then 
Xi + i = —gi — ■ ■ ■ — <72 + x\ + 1. Thus, up to this last term x\ + 1 the duality functional H(x, y) can be 
written in a similar form as G. We are not aware of a way to derive, for instance, the H duality from 
the G duality (the proof of the latter is considerably simpler). 

From the self-dualities of Theorem 12.201 we can derive four more self-dualities. 


Theorem 2.22. The J = 1 higher spin zero range process with transition operator B a ' qa is dual to the 
space reversed process with generator B a,qa with respect to G m (rh,y), G m (m,y), G„(g,n), and G„(g,n). 
Precisely, 

B a ’ qa G m = G m (B a ’ qa ) T , and B a > qa G m = G m (B a ’ qa ) T , 

where the equality holds for all matrix elements indexed by m € WG)' and y £ Y. Likewise, 

B a ' qa G n = G„(B Q ' qa ) T , and B a ’ qa G n = G n (B a ’ qa ) T , 

where the equality holds for all matrix elements indexed by g £ G and n £ WY^'. For the G dualities, 
we also require both sides to be finite. Recall from Definition \2.5\ that WG k and WYj refers to rh and 
n which are strictly ordered. 


Proof. Observe that for m,n £ Z, g £ G, and y £ Y 


[Vm\q ■ q 


[: Vm]q ■ q 


N 4 - , (y) q 

. n m— l\cfJ == — 




1-q 

■it f,-r\ _!\|t 


-NUfl _ g Nm( ^ } - q N -+ i(g} 


1-q 


4 ®The update for g according to B a,qc ‘ is from left to right, and opposite for y according to B a,qa . The rightward 
movement of the particles in g besides the k left-most does not change the value of G, and likewise for the leftward 
movement of the particles in y besides the k’ right-most. 
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\9n-l\q ' 9 


[. 9n]q ' q 


n£«) = 1 


Nt(g) = i 


N 1(9) _ q^l-l(9) 

1 -q 

-Ni(ff) g-N^G) 


i-g 


We will use the above identities as well as the dualities for G and G to prove the theorem. Since the 
proofs of all dualities are quite similar, we will only present the first one. In light of the first of the 
above identities, for m € WG^ and y E ¥ we have that G m = C G for the linear operator 


^=n 


[T]i ~ 1 

1 ~q 


Here, (T/)(m) = f(m — 1) and [T]j is the operator which acts as T in the ^-coordinate. It is important 
that we have restricted ourselves to strictly ordered rh since otherwise C would depend on the cluster 
structure of m. In light of this, the duality which we wish to demonstrate reduces to showing that as 
operators from WG k —> ¥, B a,qa C G = CG(B a,qa ) T . This follows, however, because as operators from 
WG fc —» WG k , B a ’ qa and C commute. One way to see this is to observe that C is diagonalized in 
the same basis as B a,qa . In particular, because C acts symmetrically on all variables, it follows that 
= ev(z )\H E In a similar manner as demonstrated in Corollary 12. 141 this implies the commutation. 
In fact, this deduction really requires restricting to y which are well-adapted to G. As we have seen 
earlier in the proof of Theorem 12.201 it is then possible to extend from well-adapted y to all y E ¥. □ 


3. Fusion and self duality for J € Z>i 

Fusion of /^-matrices is a representation theoretic mechanism introduced in [KR87] to construct II- 
matrices with higher horizontal spin (i.e. J E Z>i) from those with J = 1 while maintaining the 
diagonalizability of the associated transfer matrices. The procedure simplifies on Z and we provide (in 
Section [3TTD a rather simple probabilistic proof using Markov functions theory. Fusion also naturally 
provides a recursion relation in J for the higher spin i?-matrix R^ \ We record the recursion relation 
in Section 13.21 and solve it explicitly in terms of g-Racah polynomials (or regularized terminating basic 
hypergeometric series) in Section 13.31 The self-dualities proved earlier in Section 12.31 immediately 
generalize to all J E Z>i. In Remark 13.181 we observe how our R-matrix can be analytically continued 
so that q J a is replaced by an arbitrary (3 € C. We comment briefly on the implications, and develop 
this further in Section [5j 

Definition 3.1 (General J higher spin zero range and exclusion process transition operators). For 
J E Z>i defincFl 

jga,q J a jga,qajgqa,q 2 a m m # j^q J ~ 1 a.,q J a 

and likewise B a,qJa . Also define 

'-j~a,q J a 'j-'a,qa'j~qa,q 2 a # # # 'j-q J ~ 1 a,q J a 

These correspond to taking J steps of the processes with parameters a, qa ,... , q J ~ 1 a. 


19 


Though the value of ev(z) is unimportant, it can be calculated to be 


ev(z) = J] 


1 - 


1-9 
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The order does not matter since each operator is diagonalized in the same basis. 
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The left-hand side of Figure [6] illustrates the sequential composition of the B qj la ’i ja operators for 
j = 1 Since each transition operator is stochastic, their product is as well. Moreover, it follows 
from Proposition 12.121 that: 

Corollary 3.2. If \jE^F i+^ia | < 1 f or 1 < * < ^ and 1 < j < J — 1, then 

i) (g)= n 'f/f* *?(«)■ (3-i) 

l ~r OtZi 

i=l 

The of choice parameters a,qa ,... ,q J ~ 1 a implies the telescoping of the product of eigenvalues and 
hence the simple form of the result. 

The following is an immediate corollary of the corresponding J = 1 duality results contained in 
Theorems 12.191 12.201 and 12.221 along with the fact that the operators B a,qa commute for different 
values of a (and the same fact for B a,qa and F a,qa ). Recall the duality functionals from Definition 12. 161 


Corollary 3.3. For all J G Z>i we have the following Markov dualities: 

• For all matrix elements indexed by x € X well-adapted to H and all y G Y k , 

• For all matrix elements indexed by g G G k and y G Y, or by g G G and y G Y k , 

_ G [B a,qJa ) T . 

• For all matrix elements indexed by g G G and y G Y, provided both sides of the equation below 
are finite, 

ga,q J aQ _ a ^ T . 

• For all matrix elements indexed bB m G WGj and y G Y (for the G duality, we also require 
both sides to be finite), 

B a ’ qJa G m = G m (B a ’ qJa ) T , and B a ’ qJa G m = G m (B a ’ qJa ) T . 

• For all matrix elements indexed by g G G and ft G WYf (for the G duality, we also require both 
sides to be finite), 


B a ’ qJa G n = G n (B Q ’ qJa ) T , and B a ’ qJ ° G n = G n (B a ’ qJa ) 

3.1. Fusion. We can describe the update procedure for the Markov chain with the transition operator 
B a,q a started in a state g G G k in the following manner. 


Jcv\ T 


Definition 3.4. For xGZ, define the product of i?-matrices 

so that 


[(4Y 


has matrix elements 


] x ,i ■ [ V/ ]- ® P 1 ]! ® ® P'jj ->• (V T ] X ® [H 1 ]! ® ■ ■ ■ 


[{ R a ] )^ qJ ] x ,i{9x,h x ,i,--- ,h x ,J-,g' x ,h x+ gi,■■■ ,h x+ i,j) 


® [M 1 ] 


with g x ,g' x € [V^x and h X)V G [IHI 1 ]^, for 1 < y < J. These matrix elements represent the transition 
probabilities from inputs g x , h X} i, ..., h x _j to outputs g' xl h x+ \.\, ..., h x+ \j. In terms of the right-hand 
side of Figure O these provide the transition probabilities from the arrows coming into a column (such 
as the one in grey) from bottom and left, to those leaving to the top and right. 
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Recall from Definition 12. Sl that notation WGj and WY{ refers to m and n which are strictly ordered. 
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Figure 6. Left: The trajectory (in terms of arrows) of the application of the Markov transi- 
tion operator B a,qa . B qa,q 01 , through B q a,q a (here J = 7). The bottom grey row represents 
the input to B a,qa and the higher grey row the output. Right: Due to the sequential update 
rule for each B q a,q “, instead of sampling row by row, one can also sample this trajectory 
sequentially by updating one column (such as indicated in grey) at a time from left to right. 


As in Definition 12.61 we use (i?^)® 9 ^ to update sequentially in the following manner. Find the first 
x E Z such that g x > 0. In the first sequential update step, let h XjV = 0 for 1 < y < J and randomly 
choose g' x and h x+ i >y for 1 < y < J according to the stochastic matrix (Ra^® qJ . The randomly chosen 
h x +i,y become input, along with g x +i for the next column update step, and so on sequentially increasing 
x. It is clear that the above described update from g to <f agrees with B a,q a . 

This update procedure can be recast in the following notation (see the right-hand side of Figured]). 
First, for x E Z define S to be the set of basis elements of [H 1 ]! 0 • • • <g> [1 H 1 ]j and S' to be the set of 
basis elements of [in' 7 ]!. In other words, S is isomorphic to {0,1} J and S' to {0,1,... , J}. For all iGZ, 
define a Markov transition operator P x : S — >• S whose matrix elements indexed by (h Xj i,..., h X) j) E S 
and (h x+ 1 , 1 ,..., h x+ i,j) E S are given by 


Px{h x ,ii'' m ih x> j',h x - (-ip, • ,/i x +i, j) — [(-^q )) ] x,i ^3xi h x ,i: i h Xj j', g x , h x +i t i, ■ ■ ■ (3.2) 

where we consider g E G k fixed, and have set 


3x — 9x T h x h x - )-i, 

As a corollary of the above definitions we have 


J 

with h x = Y h x , y 
y =i 


(3.3) 


Corollary 3.5. Fix g € G k and assume without loss of generality that x = 0 is the smallest x € Z such 
that g x > 0. Consider the Markov chain with state space S, transition operator P x , and initial state 
ho : y = 0 for 1 < y < J. Denote the value at dime’ x as {h x , y }i< y <j■ Then, from the trajectory of this 
Markov chain we can deterministically compute g 1 via \S. S\) and the probability of having an output if 
given g is equal to B a,q a (g,if). 

Define the function (j) : S —>• S' which takes (h x> i,..., h Xi j) E S to h x = f2y = \ h x>y E S'. Define the 
operator : S —>• S' which acts on functions as {Qf)(h x ^ i, • • • , h Xi j) = f(${h x , i, • • • , h Xi j . In other 
words, the matrix elements of <h, indexed by (h x> i,..., h Xy j) € S and h x E S', are given by 

^(^,1)' ' ' > h X ,J] h x ) = ^-(j)(h Xt \,...,h Xi j)=h x • (3-4) 
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Define an operator A : S' —>• S whose matrix elements indexed by h x E S' and {h x ,\ , ..., h x ,j) E S 
are given by 

(hx,li ■ ■ ■ t — Z ^■<j>(h Xt i,...,h Xt j)=h x ' || > (3-5) 

V • h'X,y = 1 

where Z equals the sum of the weights 1 <j>(h x , u ...,h Xt j)=h x ' Yly.h x , v =i9 V over a11 (h x ,i, ■ ■ ■ ,K,j) E 5. 
Observe that A does not depend on x, and, moreover, it is a Markov transition operator, meaning that 
for each h x E S', A (h x ;-) is a probability measure in the second slot (over the set S), and for each 
(h Xt i,..., h Xt j) E S, A^-; (h x> i,..., is bounded and measurable in the first slot (over the set S'). 

Proposition 3.6. The following two identities hold: 

(1) A<f> = I, the identity operator on S', 

(2) For each x E Z the Markov operator 

Q x = A P x d> (3.6) 

from S' to S' satisfies A P x = Q x A. 

Proof. The first identity amounts to the claim that for any h x , the probability measure A (h x ;-) is 
supported entirely upon <)> _1 (h x ), the pre-image of h x under cf. This, however, is immediate from the 
definition of A. 

The second identity relies on more involved properties of the /2-matrix weights. These are recorded 
in the following lemma (see also Figure [7] for a pictorial representation of the three identities of the 
lemma). 


Lemma 3.7. For all a and all m E Z>o 

qR^\m, 0; m — 1,1 )R^2 ( m — 1,0; m — 1,0) = R^\m, 0; m, 0 )R$ ( m - 0; rn ~ 1) 1), (3-7) 

qR^\m, 1; m, 1) R'qa ( m i 1; m + 1,0) = R^\m, 1; m + 1, 0)R^(m + 1,1; m + 1,1), (3.8) 

q^R^\m, 1; m, l)R^(m,0;m,0) + qR,2\m, 0; m — 1 , 1 )Rqa(m~ 1,1; m,0)) 

= R^\m, 1; m + 1, 0)Rga (m + 1,0; m, 1) + qR^\m, 0; m, ®)Rqa (m, 1; m, 1). (3.9) 

Proof. Each identity is readily checked by direct calculation. □ 


We will use the lemma to show that for any j € {1,..., J — 1}, 

7 l'x (hxi h x +i,i , • • •, h x +i,j ; h x +i,j+ii • • • Wy) (3-10) 

— g x+i,j+i x+i,j \p ^ ^h x ] h x ~\~ i , i^ • • • j h x +\ i h x -\-i,j i j) • 

Notice that the terms h x +1 j and h x - i-ij+i have been switched between the left-hand and right-hand 
sides. 

Lemma 13 . 71 readily implies (13.101) . To see this, first write 

LHSpTTOl) = E A( h x 5 h Xi \ 5 ••• 5 Ar,j) P r (Ar,i) • • • 5 hx,j 5 ^ x -\-i,i'> • • • ? hx+i, j) • 

hx,l ,... Tx, ^ (fix') 


Let (j, j + 1) denote the permutation which transposes j and j + 1, and Id denote the identity permu¬ 
tation. Then it follows by double counting that we can continue the line of equalities as 

= 1 E 

,h x 


j 


^ ^ A(/ix! ^x,cr(l); 
CTe{ld,(jj+l)} 


> K,a(J))Px(h x , <x(l), ■ ■ • j h X ,<j{ J) i ^x+l,l> ■ ■ ■ ) ^x+1, j) • 


We may now expand the definition of P x into the product of //-matrices and factor out all terms 
unaffected by the permutation cr (i.e., terms not involving level j or j + 1). To facilitate this expansion, 
let us temporarily introduce the notation that for any permutation a, gf ^ = g x (which is given) and 
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Figure 7. Pictorial representation of the three identities in Lemma 13.71 The weight of a 
diagram is given by the product of i?-matrices for each vertex. The value of the a parameter 
is specified as either a or qa depending on the height of the vertex. 


9x,y+i = 9x,y + h X) a(y) — h x+ \ )V . Note that for a = (j, j + 1), the only value of g a (x , y) which may differ 
from the case a = Id is for y = j + 1. Then, we can continue the above line of equalities as 


2 ^ ^ \{h x■ h X) \i 

hx,l 


J 

■ ■ j h Xi j) n Ra )( y9x,y,h x ,y 1 9x,y+ 1) h X -\-\^y) 
y =i 

v¥ z hj+ 1 


x 'y 1 9 x,a 9i+ 1 ) x ’j +1 R( Ji \g x j,h x , (7 (jy,g X j +1 ,h x +i i j)R( l \g X j_i_i,h xcr (j_i_iy,g x j-\-2ih x +ij+i). 

o-e{id,(ij+i)} 


We have used here the fact that A(/t x ; h X (T yiy ..., hx, a (J)) = q hx ’ rT ^ +1) hxJ+1 \(h x ;h X! i,... ,h Xi j). 

Finally, observe that by applying Lemma 13.71 we have 

y ^ Q + ’ J+ Ra ^ (9x ,j j h X(7 ^j ^, 5xj+l) h x +l,j)Ra ^ (SnjJ+l) h Xj(7 (j-\-l) j gx,j+2y ^-x+lj+l) 

{ld,(ij+l)} 

q h x+1J+1 -h x+1J J2 q hx ^+^- h ^ +1 R^(g xJ , h X:<7(j ygl j+1 , h x+ljj+1 )R^ {g^ j+1 , h x ^ j+1 y g X)j+2l h x+1J ). 

CTe{id,(jj+i)} 


crG 


Notice that besides the factor of q hx+ the change on the right-hand side is that h x +] j and 
h x+ \ t j + i have been switched. Plugging this equality into the above line of equalities, and gathering 
terms back into their original form we arrive at 

q j:+1,j + 1 x+1 ’i KP x (h x \ /i-2+1,1, ■ ■ ■ , + h x +lj ,... h x + i t j) , 

as desired to prove (13.101) . 

It remains to use (13.101) to prove the second identity of Proposition 13.61 Observe that 


RPx {h X: h x - 1-1,1, • • • , /laj+1, j) 


y ] h.Px{h'Xi4 > {hoc+l )) j) 

hx+i €.S' 


A Px{h x 


x\ hX-\-l,ll • • • ? ^x-\-l ,J 


hr. 


h.Px(hx\$ ^(^x+i)) 


(3.11) 

(3.12) 


22 

















































We claim that 


5 


Pr 


rS 


P, 


X+l. 


■5 


P 


x-\-2 


$ 

▲ 

A 

$ 

A 

A 

$ 

A 

A 

▼ 

u/ 

T 


▼ 

_ 


Q a 


Qx +1 Qx+2 


Figure 8. The intertwining relation for Markov function theory. 


hP x (h x ', /lx+1,1, • • • , /lx+l,j) , \ /q 1 Q\ 

^ +i)) — A^x+ij Rx+i,ii • • ■: h x +i y j) • (3.13) 

This follows from two facts. First, given h x + 1 , the left-hand side of ()3.13l) is a probability measure on the 
set 4>~ 1 (h x+ 1 ) C S. To state the second fact, let us introduce short-hand that LHS (|3.13p ( JJ+1 ) is equal 
to the left-hand side of (|3.13l) when h x +ij and h x+ ij+i are switched. Likewise define RHS (13.1311 Ld+ 1 ). 
It follows from (13.1011 that for any 1 < j < J — 1 , 

LHS (l3T3l) RHS (l3T3l) 

LHS imi Od+i) - RHS IXim Oj+O' 

In fact, both sides are either equal to g , 1 or q. Thus, combining these two facts (along with the fact 
that transpositions (j.j + 1), 1 < j < J — 1 generate the symmetric group S(J)) proves (13.1311 . 

To complete the proof of the second identity of Proposition 13.61 combine (13.1111 and (13.131) to find that 

APx {^h X i (^x+iji) • • • j ^ ' AP x (h XJ <j) (Ah-i))A yh x +i: (hx+i,i > • • • ?^x+i, j)^■ 

Aa+ie5' 

But A P x (h x ; (t>~ l {h x+ 1 )) = AP x $(h x , h x+1 ) = Q x (h x ; /i s +i) implying APj, = Q X A. □ 

We make use of the following sufficient condition for when a function of a Markov chain is itself 
Markov. This comes from |RP8lj . though since we deal with an inhomogeneous Markov chain, we 
provide the (essentially unchanged) proof. 


Proposition 3.8. Consider measurable spaces S and S' and a measurable transform f> : S —l S'. Let 
$ be a Markov transition operator from S —>• S' defined according to T/ = / o <f. Consider a collection 
of (x E Z>o )-indexed Markov transition operators P x : 5 —l 5. Suppose that there exists a Markov 
transition operator A : S' S such that (see Figure m 

• AT = I, the identity operator on S', 

• For each x E Z> 0 the Markov operator Q x = A P x & from S' to S' satisfies A P x = Q x A. 

Let X(x) be Markov with x-indexed Markov transition operators P x (so that P x takes one from the x to 
x + 1 state) and initial distribution A (y, •) where y E S'. Then Y(x) = 0(X(x)) is Markov with starting 
state Y(0) = y and x-indexed Markov transition operators Q x . 

Proof. For Borel functions / : S —> R and f : S' —>• M the first condition implies thatEl A(T/')/ = f'Af. 
Using the second condition we find that 

A P x ($f')f = Q x A($f')g = Q x f'Af. 

In the same manner, we find that for any x E Z>o and Borel functions / : S —l M and /q, ..., f' x : S' —> M, 
A(T/')P 0 (<h/0 • • • P*-!($/*)/ = foQofiQiti • • • Qx-i/'A/. 

This immediately implies the conclusion of the theorem. □ 

“^Concatenation of operators and functions should be read from right to left, unless indicated by parentheses. 
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Left to Right Update 
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Figure 9. Illustration of the J higher spin zero range process (on top) and exclusion process 
(on bottom). Here J is at least three (as there are clusters / jumps of size at least three. The 
17-matrix R^' 1 determines the distribution of the sequential updates. 


We return now to the specific definitions of S, S', P x , 4>, T, A, and Q x given earlier in this section. 

Corollary 3.9. Initialize the Markov chain with state space S and transition operator P x to initial state 
ho ; y = 0 for 1 < y < J and denote the value at ‘time’ x as { h x ,y}i<y<j■ Then h x = ..., h Xt j) is 

Markov with starting state ho = 0 and Markov transition operator Q x . 

Definition 3.10 (General J 17-matrix). The J higher spin 17-matrix R J a : V 7 ®!!' 7 — > V 7 ®!! 7 is defined 
such that (here x € Z is an arbitrary and inconsequential index) 

[4 J) ki(s..W„ hx+l) — Ig x +h x =g' x +h x +iQx(,h'Xi hx+l^i (3-14) 

where Q x is defined with respect to g x (via P x ) in (13.61) . 

We may now define the higher spin versions of the zero range and exclusion processes introduced in 
Definitions 12.61 and 12.91 respectively. Observe that these depend now on four parameters, q,v,a,J. As 
before, we assume that q, n,a satisfy the conditions of (12.21) . 

Definition 3.11 (General J higher spin zero range process and exclusion process). For J E Z>i define 
the J higher spin zero range process according to the sequential construction of Definition 12.61 with 
the 17-matrix r! ( > ' > defined in (13.141) . Now, anywhere between 0 and J particles can move right in each 
update (i.e., the h x are basis elements of H' 7 ). We still write the corresponding Markov chain at g(t). 
Likewise, define the space-reversed zero range process y[t) as in the J = 1 case. Similarly, define the 
J higher spin exclusion process as in Definition 12.91 with the 17-matrix Ra ■ Now, particles can move 
right by up to and including J spaces. See Figure [9] for illustrations of these processes. 

Corollary 3.12. Recall Definition \3. 11[ The transition operator for the J higher spin zero range process 
is B a ’ q a , it space-reversal is B a,q a , and the exclusion process is p a ’ q “. 

Proof. These results follow from Corollaries 13.51 and 13.91 □ 
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Figure 10. The decomposition of a higher horizontal spin vertex. 


3.2. Recursion relation for R-matrix weights. In ()3. 14jl we have provided a method to calculate 
rP in terms of Rp . We can also write a compact recursive formula (in J ) that Rp must satisfy. 

In defining Rp we applied A to h x E so as to yield a probability measure on the pre-image 

1 (h-a;). Let p < ' J \h x .\ \h x ) denote the projection of this probability measure onto h X} \, the first outgoing 
arrow. In other words, p^ J \0\h x ) and p( J )(l|h x ) (respectively) represent the probability that h X) i = 0 
or 1 given the value of h x . One calculates that for j E {0,1}, 

^ (j\hx) = 'y ^ A (h x ‘, (j = h Xt i, h x ,2■,■■■, = 

{j = hx,ijhx,2j”’jh x j')^.(j) ^ (hx') 

Lemma 3.13. The R-matrix R.P satisfies the following recursion relation (let i 2 = i\ + j\ — j%) 
Ra\h,ji,i2,h)= P2 P ( ' J \ a \h)Ra\h,a;ii+a-b,b)R[ J a ~ 1 \i 1 +a-b,ji-a-,i 2 ,j 2 -b), (3.15) 

a,6g{0,l} 

which, along with its value at J = 1, uniquely characterizes its value for J E Z>i. 

Proof. This follows from the definitions of rP and p^ J \a\j\), and the decomposition in Figure [TUI □ 


K 1 -q 


3.3. g-Racah form of R-matrix weights. We explicitly solve the recurrence relation in Lemma 13.131 
in a form similar to that which was (to our knowledge) discovered in [Manl4| . We provide an explicit 
proof of our R-matrix formula. 

Definition 3.14. The regularized terminating basic hypergeometric series is given b\Fl 


r+lf^r 



n / r 

zk ( ■ ’ \ k II( a * ; q)n-k 

k =0 vTQh i=1 

= II(^ «)" ' r +^ r ( Q b[, , 1 .’. , br 

i= 1 v - 
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Pochhammer symbols (a; q) n , where n is allowed to be negative are defined as in Mathematica by: 


(n, q)n 


\mZp-aq k ), n> 0, 
< 1, n = 0, 

lrik"o' 1 (l-aq n+k ), n< 0. 
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Theorem 3.15. Fix J G Z>i. Then, with (3 := aq J , 


( , 2j 1 -j 1 2 j 2 -j 2 | h+h | 

R a in,Jl-,l2,32) = lu+ji=i 2 +j 2 <? 4 4 ' 4 ^ 2 


yji * 2 qJ 2 jl+* 2 ^_ Q ,j / , 


4^3 


-1 


q 12 ;q n ,—/3,—qva 

v , q 1 +-? 2- * 1 ; ( 0a _1 g 1_ * 2_ '? 2 


9,9 


(g;g)i 2 (-a;g)* 2 +j 2 (^« V jl \q)ji-j a 

In particular, R^P is stochastic and satisfies the recursion \3.15\) . 

Remark 3.16. To relate to [Bor hi, we may rewrite 

RvP(h,ji;i2,j2) = ^ ^ j q i _ jl ,q l 00 u’i J \i2,j2,iiJ l), 

(9V ^;q)oo 

where the parameters (s , v) are related to our parameters (a, v) as a = —sv, u = s 2 and v = 
—a/y/u, s = y/v. The general J vertex weight Wy from |Borl4l Corollary 6.5] is given by: 

wi J \h,ji\i2,h) 


L u+ii=i 2 +i2' 


• 2 -2 

i • x '\ ~^ x ‘2 . *1 (ll —l)+io?o ... 


ll v n (vs 1 ;q)j 1 -i 2 


(' q-,q)ii(vs-,q) il+jl 


403 


q 11 ;q l2 ,q J sv,qsv 1 
s 2 ^ g!+ii —*2 ^ ^1+ J- u — j 1 


9,9 


Proof. In order to prove this, it suffices to check that the J = 1 formula in the statement of the theorem 
matches the formulas given in Definition 12.11 and then to check that for J G Z>i, the formula in the 
theorem satisfies the recursion (13.151) . The J = 1 case is easily checked from definitions so we proceed 
to the recursion. To achieve this aim, we utilize an identity (|3.19D involving 403, as well as a recasting 
of our R-matrix in terms of q-Racah polynomials and an associated three term recursion for those 
polynomials. 

Let us start by rewriting Ra and the desired recursion. Changing from a , v to s, v variables via 
a = —sv, v = s 2 , we can rewrite 


R { a J) (h,ji,i2,32) = li 1+jl=i2+j fi-l ) i2 +i 2 +n q — 

s jl+j2-l2 y j2-jl+l2 ( VS ~ 1 ;q)j 2 _ il 


2Jl —j 2 2J2~j 2 , * 2 +if I h.O?,- 1 )+ i l h 


(3.16) 


-403 


q 12 ;q n ,q J sv,qsv 1 

s 2 , —*! q 1 + J - i 2 — j 2 


9,9 


(95 q)i 2 { vs i 9)i2+i 2 (9 J 9 1 jl ',q)ji-j 2 
For the remainder of the proof we will replace (*i, ji, * 2 ,^ 2 ) by (g. h, g ', h'). By utilizing the formula for 
Ra and p^ J \ the desired recursion can be rewritten as follows (by agreement, g' = g + h — h!\ note that 
this quantity does not change throughout the identity): 

RP\g,h,g',h') 

= _ i, 9 >' _ I, + ( !,-^, (1 ^ ) ^b + i./-, 8 '.y) 

(1 — su)(l — q J ) H (1 — su)(l — q J ) H 

sv(l - q 9 )(q J - q h ) ( j_ 1} , , , ,, _ n _ (1 - svq 9 )(q J - q h ) {J _ X) , 

+ (1 — sv)(l — q J ) ^ qa 5 *-,h,g ,h 1) (1 — sv)(l — g J ) ^ qa ^ 9 -,h,g ,h). (3.17) 


Using (I3.16p . we can recast this desired recursion as an identity between 403 functions: 
q~ 9 , q J sv , qsv~ 


403 


1 s' - a 3 a 9 sv. asv 1 \ {q h — l) (sq 9 — v) - (q 9 '-, q 9 ,q 9 sv,sv 1 


^ 1 —9 ^ q ^-~\- J — h —9 


Q,Q = 


-403 


{q J - 1) (s - v ) 

V ( q h - 1) (s 2 q 9 - 1) (q W - q J ^j q 9 + h ~ h ' 
(q J — 1) (s — v ) ( svq 9+h — 1) 


s 2 ,q h ' 9 ,q 1+J h 9 


9,9 


403 


9 g ';q 9 1 ,q J sv,sv 1 

s 2 q h '- gq J - h -g 


9,9 
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+ 


+ 


403 


q 9 (q 9 — 1) ^ sq 9 — vq h '^j 
(q J -l)(s-v) 
qh~h' ^ q h ' — q J ^j ( svq 9 — 1) ^ vq h ' — sq 9 ^j 
(q J — 1 ) (s — v ) ( svq9 +h — 1) 


q 9 '\q l 9 ,q J sv,sv 1 
s 2 , q 1+h '~ 9 , q 1+J ~ h ~9 


q,q 


403 


q 9 '\q 9 ,q J sv,sv 1 

s 2q\+h'-g q J-h-g 


q,q 


For better notation, let us set q 9 = G, q h = H , and replace q h ' by q 9 ' GH. We arrive at the following 
identity: 


403 


q 9 ';1 /G,q J sv,qsv 1 \ (1 — H)(Gs — v) 

q,q)=-n - 777 -r403 


s 2 ,Hq 1 ~ 9 ' ,q 1+J /(GH) 


q 9 ’; 1 /G, q J sv, sv 


-l 


(1 — q J )(s — v) \s 2 ,Hq- 9 \q 1+J /(GH) 


q,q 


v(l-H)(GH-q 9 '+ J )(l-s 2 G) - (q-9'-,q- 1 /G,q J sv,sv- 1 


(1 — q J )( 1 — GHsv)(s — v) 


403 


s 2 ,Hq~ 9 ,q J /(GH) 


q,q 


+ 


+ 


(1 — G)(s — Hvq 9 ) - / q 9 '■, q/G, q J sv, sv 


(1 -q J )(s-v) 


‘403 


,-l 


s 2 ,Hq 1 ~ 9 ' ,q l+J /(GH) 


q,q 


(GH - q 9 ' +J )(s - Hvq 9 ')(1 — Gsv) - fq 9 ';1 /G,q J sv,sv 


(s — i>)(l — q J )( 1 — GHsv ) 


4 03 


„-i 


s 2 ,Hq l ~ 9 ' ,q J /(GH) 


q,q • (3.18) 


Note that for integer g' = 0,1,2,..both sides of the above identity are rational functions in G, H, q J , s, v, 
and q. 

We will now use the following simple identity for basic hyper geometric series (which can be readily 
verified by looking at individual terms in 4 ^ 3 ): 


403 


q,z 


q n ] b, c, d 
U, V, w 

(1 — d)(c — w) 
(c — d) (1 — wq n ) 


4 03 


q n -b,c,dq 
u, v, wq 


, , (1 -c)(w-d) - 

< '’ z, + (c-<i)(l-u . 9 ") 403 


q n ;b,cq,d 
u, v, wq 


q,z 


(3.19) 


Using this identity, one can rewrite summands in the right-hand side of (13.181) with bottom arguments 
Hq~ 9 or q J /(GH) in terms of basic hypergeometric functions with bottom arguments Hq l ~ 9 and 
q 1+J /(GH). Let us denote the three resulting functions by 


403 


-1 


q 9 '; q/G, q J sv, sv 
s 2 ,Hq 1 ~ 9 ' ,q 1+J /(GH) 


q 9 '\ l/G, q J sv, qsv 


-1 


1 ' 4 ^ 3 Vs 2 , Hq 1 - 9 ', q 1+J /(GH) 


q,q 


q,q 


& q -l := 403 


q 9 '\q 1 /G,q J sv,q 2 sv 1 
s 2 ,Hq 1 - 9 ' ,q 1+J /(GH) 


q,q 


Applying (13.191) to the first term in the right-hand side of (13.181) . we get 


,-1 


(1 —H)(Gs — v) 1 f q 9 '',l/G,q J sv,sv 
(1 -q J )(s-v) 403 \s 2 ,Hq- 9 ',q 1+J /(GH) 


q,q = 


(1 — G) (^s — Hq 9 'v^j 
(1 - q J ) (s - v) 


^ 1 - GHq~ 9 ^ 

® q +— — y- $ 1 . 

1 — q J 


The second term in the right-hand side of (13.181) is rewritten as follows (note that we need to use (|3.191) 
twice): 


v(\ - H)(GH - q 9 ' +J )(l - s 2 G) - (q-v'-q^/G,^ sv, 


(1 — q J )( 1 — GHsv)(s — v) 


403 


, sv 


-1 


s 2 ,Hq~ 9 ' ,q J /(GH) 


q,q 
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+ 


(G - 1 )v(Gq - 1) (Gs 2 - 1) q- g ' (sq 9 ’ - Hv) (GHs - vq J ) 
( q J — 1) (s — v){Gs — v){GHsv — 1 ){Gqs — v) 

Gv{Gq - 1) (Gs 2 - 1) q~ 9 ' ( q J - H) (Hv - sq 9 ') 

( q J — 1) (Gs — v){GHsv — l)(Ggs — v) ^ 

v(Gq - 1) (Gs 2 - 1) q~ 9 ' ( q 9 ' - GHq) ( vq J - GHqs) 
( q J — 1) (GHsv — 1 )(Gqs — v) {Gq 2 s — v) 

Gv (Gs 2 - 1) q- g ' ( q J+1 - H) (qs - v ) (q 9 ' - GHq^j 




$1 


(■ q J — 1) (GHsv — 1 ){Gqs — v ) ( Gq 2 s — v ) 


V 1 - 


The third term already contains Tg, so there is no need to apply (13.191) to it. The fourth term in the 
right-hand side of (13.181) takes the form: 

{GH - q 9 ' +J )(s - Hvq~ 9 ')(l - Gsv) - {q~ 9 '; 1/G, q J sv, sv 


(s — n)(l — q J )( 1 — GHsv ) 


403 


,-i 


s' 2 , Hq 1-9 ',q J /{GH) 


q,Q 


(G - 1 )q~ 9 '{Gsv - 1) isq 9 ' - Hv) (■ vq J - GHs) Gq~ 9 '{Gsv - 1) (, q J - H) (sq 9 ' - Hv) 
V J — $> q + - V 7 


( q J — 1) (s — v)(Gs — v){GHsv — 1) 


{q J - l){Gs-v){GHsv - 1) 


d*i. 


We see that the right-hand side of (13.181) can be written as a linear combination of <hg, d>i, and $ ? -i. 
Note that the left-hand side of ()3.18|) is a multiple of <3?i. Moreover, collecting coefficients by & q , $i, 
and d'g-i, we see that (|3.18|) is equivalent to 


GHq 9 ' ^1 — q 9 ') (l — svq J ) 


Ti = B$ q -1 - {B + D )<£>i + D® q , 


(1 — q J ) (1 — GHsv ) 
where we are using the notation 

Gv {Gs 2 - 1) (q J+1 - H) {qs - v) (l - GHq 1 - 9 ') 


(3.20) 


B :=- 


D := 


{q J — 1) {GHsv — l)(Ggs — v) {Gq 2 s — v) 
(1 — G) ^s — Hq~ 9 'v) {GHs — vq J ) G{q — sv) 
(1 — q J ) (Gs — v){GHsv — l)(Ggs — v) 


To complete the proof, we aim to match this desired identity (13.201) to the following (/-difference 
relation for the (/-Racah orthogonal polynomials |KS96l (3.2.6)]: 


q n (l - q n ){ 1 - abq n+1 )R n {x) = B{x)R n {x + 1) - {B{x) + D{x))R n {x) + D{x)R n {x - 1), 
where 

(1 — aq x+l ){\ — bd(/ a;+1 )(l — cq x+1 ){\ — cdg x+1 ) 


(3.21) 


and 


B{x) = 
D{x) = 

R„(x) = 403 


(1 — cd(/ 2x+1 )(l — cd q 2x + 2 ) 
q{ 1 — ^‘ T )(l — d(/ x )(b — cq x ){ a — cd q x ) 
(1 — cdg 2x )(l — cd q 2x + l ) 


q n ,ab q n+1 ,q x ,cd q x+1 
a q, bdg, cq 


q,q 


n = 0,1, 2 


,...,N, 


(3.22) 
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with either of the following three truncations holding: 



a q; 

< bdg; 
cq. 


(3.23) 


The above quantities R n (x), x = 0, 1,2,..., N, are polynomials in q~ x + cd q x+l (of degree n), and are 
orthogonal on {0,1,2,... ,N} with a certain weight. 

In fact, for the identity (13.21[) the truncations (|3.23l) are not relevant since (|3.21l) is an identity 
between rational functions in a,b, c, d. We will utilize (I3.2ip with this understanding. 

Note also that since the bottom arguments of the hypergeometric function for R n (:r) do not depend 
on x or n, the same identity (13.211) also holds if one replaces 4^3 by 4^3 in (13.221) . 

Recalling our notation G = q g and H = q h and settin sF^ 

n = gx = g- a = s 2 /g; b = vq J ~ 9 '/s\ c = q h ~ 9 '] d = sq 9 '~ 9 ~ h /v, (3.24) 


one readily sees that 


$1 = 4^3 


q 9 ; 1/G, q J sv, qsv 1 
s 2 ,Hq 1 - 9 ',q 1 +J /(GH) 


q,q) = const • R n (x), 


= const • R re (x — 1), 4> 9 -i = const • R „(x + 1), 

with the overall constant 


const = (s 2 ; q) g '(Hq 1 9 ' \ q) g/ (q 1 +J /(GH)- q) g ,. 

The coefficients in (|3.20l) and (13.21 j) also match: 

(1 - q J ) (1 - GHsv)B = GH ■ B(x)- (l - q J ) (1 - GHsv)D = GH ■ D{x). 

This implies the desired recursion on the R matrix. □ 


Remark 3.17. From the above proof we see that the vertex weights R^'P (i\. j \; * 2 , J 2 ) can be interpreted 
as an analytic continuation of g-Racah polynomials. Namely, if v = s 2 = q~ I for some I = 1,2,..., 
then Ra\iii 31 G 2 G 2 ) is a multiple of the g-Racah polynomial Rj 2 (*i) ( 13 . 221 ) on i\ € { 0 , 1 ,... , /}. The 
parameters of the polynomial R * 2 (i 1 ) arcFl 

a = q~ I ~ 1 , b = -aq I+J ~ i2 , c = q jl ~ i2 , d = -a~ 1 q- I ~ ja 


(these are the same parameters as in (|3.24l) expressed through a = — sv and v = s 2 = g -/ ). The fact 
that ag is a negative integer power of q ensures that the index *2 and the variable i\ of the g-Racah 
polynomials are restricted to the finite integer segment {0,1 

Note that the parameters of the g-Racah polynomials also depend on (ii, ji; *2 > 42 )) and one cannot 
interpret the weights as values of a single g-Racah polynomial R„(x) (with constant parameters a, b, c, d) 
at various integer points x. This observation agrees with the fact that we can choose vertex weights 
to be non-negative (see Proposition 12.31) . but the orthogonal polynomial R n (x) must change sign as a 
function of x. 


Remark 3.18. Observe that our R-matrix formula in Theorem 13.151 admits an analytic continuation 
in the parameter f3 := aq J . For f3 € C which is not equal aq J for some J € Z>i, it is unclear whether 
Ra is stochastic. One can argue, as in the comments after [Borl4t Corollary 6 . 6 ], that for fixed i ], j \, 
the sum over i 2 ,32 of the R-matrix elements equals 1 (though the nature of this convergence may not be 
so straightforward). It is the positivity, however, which becomes suspect for general /3. By inspection 

24 In fact, it is also possible to assign other values to a, b, c, d which match all the parameters. We will not discuss all 
possible ways of choosing a, b, c, d. 

"^Similarly one can take R ;i (*2), which involves interchanging parameters (a, b) with (c, d). 
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of terms in the summation defining 4(^3 one finds that unless j3 takes the special form (which results in 
many terms being zero in the sum), the terms in the summation are not all of the same sign. Empirical 
computer testing indicates that for generic /?, R-matrix entries are not always positive. 

It is because of this lack of stochasticity that we do not pursue this j3 continuation further at this time. 
One can certainly (as in the comment from [Bor 14] ) extend the eigenfunction relations to general /3. It 
may also be possible to extend the duality and the formulas proved in Section Q] to this non-stochastic 
setting (in which case, expectations in Section U] must be replaced by signed expectations). The only 
case in which we deal with an extension to general /3 is in Section 15.61 


4. Moment and c^-Laplace transform formulas 

Combining the first duality result of Corollary 13.31 with the integrability of B a,qJa given in Corollary 
13.21 we may employ results of [BCPS14] to compute moment formulas and then follow the general 
scheme developed in [BC141IBCS12] to arrive at an e^-Laplace transform formulas for the J higher spin 
exclusion process (recall from Definition 13.II) with step initial data. The other dualities of Corollary 13.31 
may prove similarly useful in computing moment formulas (perhaps for different choices of initial data), 
but we do not pursue these here. 


Theorem 4.1. Fix J E Z>i, f3 = aq J and consider the J higher spin exclusion process x(-) with step 
initial data Xj(0) = —i, i € Z>i. For all n\> ■ ■ ■ > > 1 and t E Z>o, 


E 


Y[q Xn iW +ni 


L i— 1 


, k(k — 1) 

(-1 ) k q ^^ 


(27ri 




Tl 


n 


%A - ZB 


n 


Ik 1 <A<B<k ^ ZB j-\ 


1 — VZj\ n i (\ + @Zj\t 


1 — Zj 


dzi 


l + azjJ Zj (l-vzj) 


where the simple closed integration contours 74 ,...,are chosen such that they all contain 1 . the 7 ^ 
contour contains q*y B for all B > A and all contours exclude 0 and l/v. 

Proof. This is proved in the same manner as |BCPS14i Corollary 5.19], replacing the eigenvalue Jl5=i 

therein by the present value rij=i 15757- Rather than repeating the proof, we briefly explain the idea 
behind it. The first duality result of Corollary 13.31 implies that 


E 


K 

[R 

i= 1 


x ni (t+l)+rij 


_ ga,q J a^ 


K 

[R 

2—1 


.Xni (t)+rii 


where B a,q a acts on the n variables via the associate of ft with y, as in Definition 12.51 The operator 
B a ’i a is diagonalized (Corollary 13.21) by the eigenfunctions introduced in Appendix [A] and utilizing 
the associated Plancherel theory, we arrive at the solution to the time evolution equation satisfied by 


E 


nti q Xn S t+1 l +ni , thus proving the theorem. 


□ 


For the exclusion process with step initial data, the formula from Theorem 14.11 provide a complete 
characterization of the distribution of x(t). This is because each random variable q x ' i -( t ) +n , 1 < n < N, 
is in (0,1] for all t. The knowledge of all joint moment thus suffices to characterize the joint distribution. 
Despite this fact, it is not obvious how to extract meaningful asymptotic distribution information from 
these formulas. In the case of one-point distributions (i.e. the distribution of x n (t) for a single n) 
this was achieved in [BC14] , We will apply the approach developed therein (in particular, the general 
restatement of the calculation in [BC14i ] which can be found in [BCS121 Section 3]). 

Theorem 14.21 provides two Fredholm determinant formulas for what is called the e g -Laplace transform 
of the observable g x »( t )+ n ; anc [ consequently for the one-point distribution of x n (t) (see [BC14 , Proposi¬ 
tion 3.1.1] or |BCS12l Proposition 7.1]). This type of Fredholm determinant formula (in particular that 
of (14,11) 1 is quite amenable to asymptotic analysis - see for instance |BC141 IBCF121IBCFV141IBCR121 
IBar141 lFVT3l IVetT4l [HOT4l K1SS14] . 
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Theorem 4.2. Fix J E Z>i, (3 = aq J and consider the J higher spin exclusion process x(-) with step 
initial data Xj(0) = —i,i € Z>i. For all n E Z>i, t E Z>o, and £ € C\R+, 


E 


(Cq Xn( - t '> +n -,q) c 


= det (I + K ( ) 


(4.1) 


where det (/ + K^j is the Fredholm determinant of : L 2 (C\) —» L 2 (C\) for Cj a positively oriented 
circle containing 1 with small enough radius so as to not contain 0, 1/q, and 1/v. The operator Kq is 
defined in terms of its integral kernel 

^ r ioo+1/2 


( r lc 

K < iw ' w ">=m L 


7r 


:(-CY 


l 


with 


g (w) = 


—ioo+l/2 sln V — TTS) 

(ow-q) 

oo \ n (q ) c 


g{w) _ 

g (q s w) q s w — w' 


ds 


(w,q) c 

The following second formula also holds: 


E 


1 


(£ g s„(t)+n ;g ) c 


(- aw;q)oo) {vw\q) ^ 

_ det (I + (K) 

(Ci q)o o 


(4.2) 


where det (/ + (K) is the Fredholm determinant of £ times the operator : L 2 ((7o,i) —» L 2 (C'o,i) /or 
6 * 0,1 a positively oriented circle containing 0 and 1 (but not \/v or —1/a). The operator K is defined 
in terms of its integral kernel 


K(w, w') 


g{w)/g(qw) 
qw' — w 


where the function g is as above. 


Proof. This type of deduction of Fredholm determinant formulas from g-moment formulas has appeared 
before, cf. [BC141IBCS121IBCF121IBC131IBCFV14] , Hence, we provide only the steps of the proof, 
without going into much detail. We also do not recall the definition of Fredholm determinants, cf. 
|BC14i Section 3.2.2], 

In order to prove the first formula (which is sometimes called a Mellin-Barnes type formula) we utilize 
the formula for E[g fc ( x ™( t )+ n )] from specializing all n* = n in Theorem l4.ll Call pk '■= E[g fc ( x "( t ) +ri )] so as 
to match it with the formula present in [BCS121 Definition 3.1], subject to defining f(w) := g(w)/g(qw), 
with g from the statement of Theorem m We may then apply |BCS12l Propositions 3.3 and 3.6] with 
the contour Ca = C*i, and D^d;k specified by setting R = 1/2 (and d arbitrary, as it does not 

matter for this choice of R ). The output of these propositions is that 

/-k 

E^jh = de t (/ + *c)- 

fc >0 q ' 


In the course of applying these propositions, it is necessary to check that a few technical conditions on 
the contours, as well as £ and g are satisfied. These are easily confirmed for |£| small enough, and C\ 
a small enough circle around 1. The only condition depending on the function g is that |g(i«)/g(i/ s 'w)| 
remain uniformly bounded as w E C*i, k E Z>i and s E Dji^-k varies. This is readily confirmed for g 
from the statement of Theorem 14.21 

Now, observe that for £ with |£| small enough, we also have that 


£ fc 


= E 

k>0 q 


(Cq x ^ t l +n -q) c 
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This is justified (as in jBC14l Theorem 3.2.11]) by the fact that q x ^ t )+ n g (0,1) and an application of 
the g-Binomial theorem. This establishes (14.11) for |£| sufficiently small and since both sides are analytic 
in C \ M + the general Q result of (14.11) follows via analytic continuation. 

The Cauchy type formula (14.21) follows from |BCS12l Proposition 3.10] along with a small amount of 
algebra. The proof essentially follows that of jBC141 Theorem 3.2.16]. □ 

Remark 4.3. Setting g t {i) := ®j_i(t) —Xj(t), g(t) evolves as the zero range process with step initial data 
corresponding to having <7i(0) = +oo and gi(0) = 0 for i > 1 . Let C s (t) = gi{t ) be the number 

of particles of g(t) strictly to the right of site s at time t. Then clearly { x n (t ) + n > s} = { C s (t ) > n} 
and hence Theorem 14.21 provides an exact formula characterizing the distribution of C s (t) as well. 

5. Generalizations, specializations and degenerations 

5.1. Case (2) of Proposition I2T31 Recall that in case (2) of Proposition 12.31 we assume q G (—1,0], 
a G (0, l/|g|), and v G ( — l/|g|, min(l, a/|g|)). As already observed in Proposition 12.31 under these 
conditions one readily confirms that the R-matrix is stochastic. Thus, the zero range and exclusion 
processes built from the R-matrix in Definitions 12.61 and 12.91 remain valid. The q -Hahn operators R a 
and R a remain well-defined. All of the duality and fusion results extend as we now explain. 

In Section [2] all of the results can be extended. Proposition 12.121 holds since the proof only relies 
upon |Borl4„ Corollary 4.5 (i)] which remains valid in such an extension of parameters. Alternatively, 
one can observe that the equality demonstrated by the proposition can be analytically continued into 
the regime of parameters we are considering. In order to show the analyticity, one uses the bound that 
| | < 1 which ensures the convergence of the left-hand side of equality in the proposition. The 

factor Y—vz- comes from the eigenfunctions whereas the factor comes from the weight associated 
with (0,1; 0,1), which is the only vertex weight which can occur an unbounded number of times in 
the application of the operator B a,qa . Proposition 12.131 holds since both sides are rational functions 
in the parameters and can be extended to the desired range. Corollary 12.141 holds since the results of 
Appendix [A] hold for this present case of parameters. Theorem 12. 191 holds since the only properties the 
proof relies upon are those notes in Remark l2.11l (both of which remain valid under the present choice of 
parameters). Theorem 12.201 holds since the proof only relies upon [Corl4l Proposition 1.2] which holds 
as rational identity in the parameters q, a and clearly extends to the range of parameters presently 
considered. Theorem 12.221 holds since the proof only relies upon Theorem 12.201 and the eigenfunction 
relations which remain valid. 

All of the result of Section [3] besides those pertaining to the exclusion process hold as long as the 
triple (g, v, aq 3 ) satisfy the conditions of case (2) of Proposition 12.31 for 1 < j < J. This is because they 
all follow from results of Section [2j 

We do not pursue modifying the results of Section [H As it stands, the formula in Theorem 14.11 is 
not well adapted to taking q negative (because of the nested structure of the contours). Instead, one 
might first shrink the nested contours to all lie around 1 (cf. [BCPS141 Proposition 3.2]) and then argue 
by analytic continuation that the formula in which the parameters have been extends to case (2) of 
Proposition 12.31 Such a formula would then serve as the input to establishing a result like Theorem 14.21 
We leave the justifications necessary to prove such formulas to future work. 

5.2. R-matrix reflection and inversion. Let /, J G Z>j. Temporarily, to emphasize the role of I 
and J, let us write 

? JD ^ 2 ? J 2 Ji 0 ) \ 

where the right-hand side implicitly depends on q as well as I via v = q^ 1 . 

There are two actions on a vertex we consider. The first is reflection in the diagonal under which 
(i\. ji] i r 2 , j'2) '-t O'l, *i; J 2 , *2) an d the second is inversion of arrows under which (*i, ji] 12, J2) l_ t(^ — 
ii, J — ji; I — 12, J — ^2). The result of these two actions on the R-matrix is quite simple. 
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Proposition 5.1. Let a = 1/a and q = 1/q then 

R(ji,h;j 2 ,i 2 \I, J,a,q) = R(h,ji;i 2 ,j 2 \J,I,a,q), 

R(I -h,J-jr,I-i 2 ,J~j 2 \I,J,a,q) = R(h,ji;i 2 ,j 2 \I, J,&, q). 

Composing the two transformations (reflection and inversion) results in yet a fourth stochastic R-matrix 
in which (I, J) i->- (J, I) and a , q remain fixed. 

Proof. These relations can be confirmed directly from our explicit formulas for the R-matrices given in 
Theorem 13.151 □ 

Cases (3) and (4) of Proposition 12.31 are related by inversion of arrows: The range q £ [0,1) maps 
to q € (1, Too) and the range a < —q^ 1 maps to — q~ I < a < 0. Despite this relationship, we will 
still consider each case separate. This is because when we construct our zero range process from the 
.R-matrix, the inversion of arrows takes a state g with ^ = k to an state with an infinite number of 

particles. 

Note that applying reflection to cases (1) and (2) of Proposition 12.31 leads to different R-matrices and 
corresponding processes than we have presently considered. We do not pursue the study of the resulting 
systems any further here. 

5.3. Case (3) of Proposition I2T31 Recall that in case (3) of Proposition 12.31 we assume that q £ (0,1), 
u = q for I £ Z>i, and a < — v. With this choice of /, Y 1 is finite dimensional and given by the span 
of {0,...,/}. Consequently we must replace G,G k , WG k , ¥, Y k , W¥ l " by G/, GJ, WGf, Yj, Y k , WYj 
as given at the end of Definition 12.51 The construction of the zero range processes with transition 
operators B a,ga and B a,qa remains unchanged. However, in this finite spin case, we cannot define the 
exclusion process with a right-most particle (whose transition operator was This is because 

that construction required having an infinite gap and in the finite spin setting, that is not allowed. It 
may be possible to define an exclusion process with a doubly infinite state space, so long as the gaps 
between particles is bounded by I. However, as we have no use for that presently, we do not pursue it. 
Definition 12.11 remains well-defined, even though (p q - a ^ v (s\y) is not necessarily positive (in fact, it has 
sign (—l) s ) for integers 0 < s < y < I. Likewise, the two equations of Remark 12. Ill remain valid (one 
readily observes that 1 + aq y ~ s remains non-zero for our choice of a. and for all 0 < s < y < I). 

Appendix [A] describes how the Plancherel theory is modified for the present choice of parameters. 
The eigenfunction relation of Proposition 12.121 follows since the results of |Borl41 Corollary 4.5 (i)] are 
stated in sufficient generality. It is also possible to show that the left-hand and right-hand sides of the 
relation in the proposition are analytic functions of v and a in a suitable open domain which connects 
the parameters from case (1) of Proposition 12.31 to the present choices. Analytic continuation then 
implies the extension of the relation to the present case. This argument only requires showing that for 
some fixed z. L , there is an open domain of (is,a) £ C 2 connecting case (1) and case (3) of Proposition [2731 
such that the relation | < 1 is preserved throughout the domain. This fact is readily checked. 

Proposition 12.131 holds as written. 

We will not make use of Corollary 12.141 (though it holds with replaced by Wf. max from Defi¬ 

nition [A3]) . We will also not make use of Definition 12. 171 and there appear not to be natural examples 
of initial data satisfying the well-adaptness of Remark 12.181 Since we have not formulated an exclusion 
process for these parameters, we do not have an analog of Theorem 12. 191 presently. Theorem l2.20l how¬ 
ever, holds for the choice of parameters. Since we cannot make use of well-adaptedness, we demonstrate 
this result by analytic continuation. Let us focus on establishing the first of the dualities (the other 
follows similarly) 

ga,qctQ _ G( i g“.9 Q ) T i 

Similarly to as above, note that there exists an open domain of (is, a) € C 2 connecting case (1) and case 
(3) of Proposition 12.31 along which < d for some 5 < 1. Consequently, both sides of the above 
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identity are analytic functions of v and a and the result follows. The reason for analyticity is due to 
the fact that for g and y fixed, the only R-matrix weight which can be used an unbounded number of 
times is corresponding with (0,1; 0, l)-vertices. That this weight is bounded above by 5 < 1 in 
magnitude ensures the convergence of the left-hand and right-hand sides and hence the analyticity. It 
should be noted that as we deform parameters, the transition operators are no-longer stochastic and it 
is only after arriving at the terminal locations for parameters that this property is restored. Theorem 
12.221 holds true by similar considerations. Notice that even if we were to try to follow the proof given 
earlier in the text, we would run into a problem if I < k. 

All of the result of Section [3] besides those pertaining to the exclusion process hold as long as the 
triple (q, v, aq 3 ) satisfy the conditions of case (3) of Proposition 12.31 for 1 < j < J. This is because they 
all follow from results of Section El 

5.4. Case 4 of Proposition 12.31 Recall that in case (3) of Proposition 12.31 we assume that q £ (1, oo), 
u = q for I € Z>i, and — v < a < 0. With this choice of /, V 1 is finite dimensional and given by the 
span of {0,... , /}. Consequently we must replace G, G k , WG fc , ¥, Y k ,WY k by G/, Gj, WG^, Yj, Y k , WY$ 
as given at the end of Definition 12.51 The construction of the zero range processes with transition 
operators B a,qa and B a,qa remains unchanged. However, as in Section r5.3l we cannot define the exclusion 
process with a right-most particle. Definition 12.11 remains well-defined, even though < p q - a ^{s\y ) is not 
necessarily positive (in fact, it has sign (—l) s ) for integers 0 < s < y < I. Likewise, the two equations 
of Remark 12.111 remain valid. 

Appendix [A] describes how the Plancherel theory is modified for the present choice of parameters. 
The eigenfunction relation of Proposition 12.121 follows since the results of [Borl4i Corollary 4.5 (i)] are 
stated in sufficient generality. Proposition 12.131 holds as written. We will note make use of Corollary 
12.141 (though it holds with W^ ax replaced by Wj. max from Definition IA.5I) . We will also not make use 
of Definition 12.171 Since we have not formulated an exclusion process for these parameters, we do not 
have an analog of Theorem 12.191 presently. 

Theorems 12.201 and 12.221 follow from analytic continuation of their analogous results in case (3) with 
the following modification: In the case of the G,G n ,G m dualities, one should restrict to g € G k and 
y £ Y k for some k,k' € Z>i (otherwise terms may fail to be finite). We may also, without loss of 
generality, assume that g € G k and y £ Y k ‘ for some k,k' € Z>i when considering the G,G n ,G m 
dualities (see the end of the proof of Theorem 12.201 for an explanation of how to go from this to general 
g and y). So, given g £ G k and y £ Y k let us see how this analytic continuation works. Let us focus 
on establishing the first duality (all others follow similarly) 

ga,qaQ _ 

This hold for case (3), in which q £ (0,1) and a < —q~ k - Think of q, a as complex variables and 
fix v = q ~ 1 for I £ Z>i. We want to extend this to also hold for case (4) in which q £ (l,oo) and 
a £ (—q~\ 0). We do this by observing that there exists an open domain of (q , a) £ C 2 connecting case 
(3) to case (4) along which \jt^\ < 5 for some 6 < 1. Consequently, both sides of the above identity 
are analytic functions of q and a and the result follows. The reason for analyticity is due to the fact 
that for g and y fixed, the only R-matrix weight which can be used an unbounded number of times is 
corresponding with (0,1; 0, l)-vertices. That this weight is bounded above by 5 < 1 in magnitude 
ensures the convergence of the left-hand and right-hand sides and hence the analyticity. 

All of the result of Section [3] besides those pertaining to the exclusion process hold as long as the 
triple (q, v, aq 3 ) satisfy the conditions of case (3) of Proposition 12.31 for 1 <3< J. This is because they 
all follow from results of Section El 

5.5. Stochastic six-vertex model. If I = J = 1 (so v = 1 /q), then the zero range process g(t) we 
have been considering degenerates to the stochastic six-vertex model. In that case, the R-matrix has six 
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non-zero vertex configurations, and can be parameterized by two parameters (b \, 62 ) which are between 
zero and one |GS92| . |BCG14| : 



Here ci = 1 — b ±, C 2 = 1 — 62 . 

The I = J = 1 vertex weights are nonnegative if and only if the parameters (a, q ) belong to one of 
the following two families: 

(1) q G (0,1), a < - 1/q ; 

(2) q G (l,oo), a G (-1/^,0). 

Case (1) above corresponds with case (3) of Proposition 12.31 whereas the case (2) above corresponds 
with case (4) of that proposition. The ratio 62/^1 = l/r? is denoted by r in [BCG14]. Under the present 
choices of parameters, Theorem 12.121 degenerates to match the eigenrelations proved in [BCG 141 Section 
3.3] via coordinate Bethe ansatz (see also [Lie67j ). In |BCG14j . the eigenrelations are used to compute 
transition probabilities for the finite particle stochastic six-vertex model. The authors then focus on the 
case when r < 1 (in other words, case ( 2 ) above where q € (l,oo)) and when g{t ) is started from the 
step initial data (<7i(0) = l;>i). Building on a combination of the approaches developed in |TW08j and 
[BGl4| . they compute a contour integral formula for E[t ln "(^)], L G Z>i, and eventually utilize this to 
compute a e r -Laplace transform formula for 

Theorems 12.201 and 12.221 yield self-dualities of the stochastic six-vertex model which, to the best of 
our knowledge, appear new. As described in IK11 . Section 2.2], as 61,62 —>• 0 (with r fixed), the 
stochastic six-vertex model has a limit to the asymmetric simple exclusion process (ASEP). |Sch97] (see 
also |IS111 BGS12] ) discovered two self-dualities for ASEP. These arise as limits of two of our dualities, 
those involving G and G n . Focusing on the G n duality, notice that 

G n (< 7 ,n) = \\g ni ■ r N "i (9) 

i= 1 

where we have used the fact that for g G {0,1}, [g\ q = g. In the ASEP limit, this duality becomes that 
recording in [BGS121 Theorem 4.1], The G dualities limits to that of I K'S 12 . Theorem 4.2], 

Dualities provide closed systems of differential (or in discrete time, difference) equations satisfied by 
certain expectations of observables. For instance, taking the expectation over the evolution of the g(t ) 
process, the duality involving G n implies that 

E[G n ( 5 (t + l),n)\ = £“’^E[G n (s(t),n)]. 



with step initial data, IK'S 12 developed this approach. We expect that a similar approach can be 
implemented to reprove the formula of |BCG14| for E[r iNT1 ^] (for step initial data) by utilizing our 
duality. The main challenge in this pursuit is to compute the direct transform T q ' v of E[G n (^(0), n)]. 
This is known for ASEP, and one may observe that the ASEP eigenfunctions are (after change of 
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variables) the same as those considered hereirPl Thus, what is known in this respect for ASEP likely 
follows for the stochastic six-vertex model. We do not pursue this any further here. 


Remark 5.2. If / = J = 2 and q,a are as in cases (3) and (4) of Proposition 12.31 then the A-matrix 
is stochastic. In this case, there can be up to two particles per site and up to two particles can move 
in each update step. If a = —1 /q 2 then the zero range process becomes a deterministic shift. In the 
six-vertex case (I = J = 1) this occurs for a = —1/q. In that case, ASEP arises from an expansion 
around this shift, as ct = — 1/q + e. We can perform the same expansion in the I = J = 2 case, setting 
a = — 1/q 2 + e. The below table is calculated from the J = 2 case of Appendix [Bj The overall matrix 
(three by three) has rows and columns indexed by j\ and 72 in { 0 , 1 , 2} (respectively) and each matrix 
entry (a length three column vector) has entries indexed by i\ in {0,1,2}. All terms of order smaller 
than e are left out of the matrix. A quick inspection reveals that subtracting the shift (the order one 
terms) leaves something which is not stochastic. Indeed, both g 3 (l + q)/(q — 1) and g 5 /( 1 — q 2 ) arises 
in this e expansion and their respective signs will always differ. 



/1-A \ 

q A — 1 


V 


/; 

q 2 — l 

0 


/ 


/ l _ 9 3 (2g+l)e \ 


q 2 -1 
g 3 (g+l)e 

9-1 

q 5 e 

1 —q 2 


1 — 


i-<r 


+ 1 


q 2 e 


{q A +q 


q 2 -1 


i -q 


g 4 e 

q 2 — 1 

0 

1--&- 


9 4 e 

q 2 — l 


+ 1 


q 2 -1 


/ 


TTT 


\ \ 


g(g+l)e 

g -1 

i q (g+2)e 

\ L - q^T 


( 


0 


1 - 


J2 c- 

\ + i 



This seems to be a negative indication as to whether one can extract higher spin versions of ASEP in 
this manner. 


5.6. g-Hahn processes. In Remark 13.181 we observed the possibility to analytically continue our R- 
matrix weights so as to depend on parameters a and /3. However, for general /? / ag J , J E Z>i, these 
weights were not always positive. The following proposition provides an exception to this, in which our 
R-matrix weights reduce to g-Hahn distribution weights (Definition 12.1011 . 

Proposition 5.3. Let a = v and f3 = —/r for 0 < v < ji < 1. Then 

(H! 3l i*2 5 J 2 ) lii+j'i=i2+J2 ’ ^juO^Ih) 

for i'l-jyi € Z>o- Consequently, the space reversed higher spin zero range process y and higher 
spin exclusion process for these parameters coincides with the q-Hahn zero range vrocess 21 \ and q-Hahn 
TASEP studied in |Povl31 ICorl4| . 

Remark 5.4. For the above choices of a,/3, the proposition implies that the i?-nratrix weights do not 
depend on j\ and hence the processes constructed from these weights become parallel update. 

Proof. This reduction is essentially proved in [Borl4| Proposition 6.7] by studying how the hypergeomet¬ 
ric functions specialized with these parameters. The equality of the associated processes constructed 
from these R-matrix weights follows immediately. Note, however, that as far as the equality of the 
processes is concerned, this could be shown directly by noting that the eigenfunctions and eigenvalues 

2 ®The limit 61,62 —^0 (with r fixed) involves scaling a, but not changing q (which is 1 /t) and since q, v are the only 

parameters which enter the eigenfunctions, the ASEP and stochastic six-vertex eigenfunctions are the same. 

"^Also called the g-Hahn Boson process, or ( q , p, i/)-Boson process. 
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coincide under this parameter specialization (cf. Proposition 12. 131 and Corollary 13.21) . This route would 
require extending Corollary 13.21 to general (3 in the manner described in Remark 13.181 □ 


5.7. Inhomogeneous parameters. It is possible to define a time-inhomogeneous versions of the zero 
range and exclusion processes considered earlier by replacing the time t transition operator (from the 
state time t to that at time t + 1) by (or likewise T As long as we assume /3t = q Jt a.t 

for Jt € Z>i, this is stochastic and hence generates a Markov chain. The duality of Corollary 13.31 
clearly extends, and since these operators for different t are still diagonalized in the same basis (owing 
to Corollary 13.21) . we are able to develop analogous results to Theorems 14.11 and 14.21 The only change 
in Theorem HU is the replacement 

/ I + /3zj \t ^ -r-r 1 + PsZj 
VI + azj) *-\l + a s Zj' 

and the corresponding change in Theorem 14.21 is the replacement in the function g 

( (~Pw,q )oo V ^ TT (-PsW,q)oo 
V (-aw, q )oo J ^ (~a s w, q) oo' 

It is also natural to consider spatial inhomogeneities. In that case, the R-matrices used to construct 
the zero range and exclusion processes should depend on location or particle number (respectively) in 
so far as that the parameter a can be replaced by a x . Under this generalization, the eigenfunctions 
considered in Appendix[A]no-longer suffice for diagonalization and no suitable replacements are presently 
known. An inspect of the duality proofs as well as the proof of the fusion procedure seem to suggest 
that these results can be modified to apply in this setting. We do not pursue this direction any further 
here. However, it is worth noting that when v = 0, these processes relate to those considered in [BC13] 
wherein duality and moment formulas were proved for both time and space inhomogeneities. The 
moment formulas were generalizable in such a manner due to the connect between the v = 0 case and 
the theory of Macdonald processes [BC13 . Section 6.2], 


Appendix A. Bethe ansatz eigenfunctions 

We recall results about the Bethe ansatz eigenfunctions, most of which come from [BCPS14] . 

Definition A.l. Assume that q,v are as in the first two cases of Proposition 12.31 Recall WY fc from 
Definition 12.51 and let W k equal the space of all compactly supported functions from WY fc to C. Define 
left and right Bethe ansatz eigenfunction 




z a(A) z a(B) ^ 


n 


1 - Mfl N 


- VZ n 


C?) 


= e n 

o-eS(fc) 1 <B<A<k 

= (-i)Ri - £ n MA) n 

•T / A \ % / T~»\ 


cr eS(k) l<B<A<k 


7(A) 


(B) ^ 


1 ~ Z, 


U) 


1 — vz r 


O') 


where z = (z\, ..., Zk) € (C \ {1, v 1 }) fc , and m q , v (n) is given by 


M (n) 


^q,u(n) = n 




J-=1 («; ' 

where ci, ..., cm(H) are the cluster sizes of n (i.e. n\ = ■ ■ ■ = n Cl > n Cl+ \ = • • • n Cl+C2 > • • • 


(A.l) 


(A.2) 


2 ®See Propositions 12.121 and 12.131 which describe certain operators for which these are eigenfunctions. 
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Let J- q,u be the direct transform which takes a function / € W k in the spatial variables n and produces 
a function in the spectral variables z according to 

(j-"7)o?)= £ nmun)- 

nG WY fc 

The function (F q ' v f)(z) is a symmetric Laurent polynomial in (1 — Zj)/( 1 — vzf), j = 1 ,k. We 
denote the space of all such Laurent polynomials by C k . 

Let J q ' v be the inverse transform which maps Laurent polynomials G € C k to functions in W fc 
according to the following nested contour integration formula: 

k 

{J q,v G){n) =/.../ dm (lfc) (z) f[ 1 --T|(n)G'(z). 

J~t - =1 l 1 Tm 1 

The contour 7 is a circle containing around 1, not containing i / —1 and such that it contains its image 
under multiplication by q. The term 

1 (—1) ( 2 * V(i ?) 2 dzj 

dm ( 1 ‘)<-> = iT Xl iH {zi-qz S ) |w’ ( ' 3) 

where V(z) = Eli <i<j<k( z i — z j) is the Vandermonde determinant. 

The following result comes from |BCPS14l Theorems 3.4 and 3.9]. 


Proposition A.2. Assume that q,v are as in the first two cases of Proposition [23 The transforms 
T q ' v and J q,v are mutual inverses in the sense that J q ' u J~ q ' u acts as the identity on W k , and T q ' v J q ' u 
as the identity on C k . 


If q , v are as in the last two cases of Proposition 12.31 some modifications are necessary. Let I € Z>i 
be such that v = q~ l and recall WY k (Definition 12.51) . the space of n = (n 1 > ■ ■ ■ > rik) with clusters 
of equal elements of length at most I (e.g. if I = 1, the equality is strict). Similarly, denote by W k the 
space of all compactly supported functions from WYj to C. When we substitute v = q^ 1 , m qM {h) (and 
hence Af?(n)) becomes supported on WY k . This is because the factor (zv; q) Cj is zero if any cluster Cj of 
n exceeds length I. Therefore, in the direct transform, the summation over WY k can be replaced by a 
summation of WYj. The inverse transform similarly requires some modificatioiF^ 


(J^G)(n) = £ 

Ah !k J y 


*w 


dm A (w) 


H (.w j ;q)\ j (vw j ;q) x . 


■^woX^Giw o A). 


Some notation used above must be explained. For a oartitiorFl A \~j k if A h k and each A j < /. Define 

(A.4) 


j , ^ (l-^) fc (-l) fc 9 2 j , 

dm\{w) := -;-;-det 


rrzi !r?T,2!... 


-1 -1 e(x) tW a 2 , 

1 TT A,' fldWj 

[[wdq2 


i w iq Xi — w j J ij=i 1 


2vri’ 


where w = (w \,..., w^x)) € and for A h k, we use the notation 


wo \-= ( Wl ,qwi,...,q Xl 1 w 1 , w 2 , qw 2 , ■ ■ ■ , q M 1 w 2 . ,wx tw ,qw Xew , ■ ■ ■ ,q X * w 1 w Xe(x) ) € C 


V(A), 


(A.5) 


The contour <■/ is a small circle around 1 which does not contain v 1 and does not contain its image 
under multiplication by q. 


2 ®When I Z, such a form of the inverse transform is also valid, with h/ simply replaced by h. 

^That is, A = (Ai > A2 > ... > 0), A; £ Z, with Ai + A2 + ... = k if A b k. The number of nonzero components in A 
will be denoted by l{\), and rrij will denote the number of components of A equal to j so that A = i m i2 m2 .... 
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Proposition A.3. Assume that q, v are as in the last two cases of Proposition \2.A The transforms 
F q,v and J q,u are mutual inverses in the sense that J q ' v F q ' u acts as the identity on W k , and 
as the identity on C k . 

Proof. The case when q € (0,1) follows by applying the arguments of [ B CP Si 4 , Theorems 7.2 and 7.3]. 
One then observes that everything is analytic in q away from 1 and thus the q > 1 result thus follows 
readily. □ 

It is useful to extend the space W k to include non-compactly supported functions which still have 
nice growth properties. In particular, for c,C> 0 define / c c ^ as those functions / : WY fc — > C 

such that |/(n)| < Cexplc^*^ n;} for all n € WY fc . Likewise define W ^. exp / cC y 

Proposition A.4. There exist c > 0 small enough and C > 0 large enough such that J q ' v F q ' u acts as 
the identity on (and likewise for W/. exp ( c j- 

Proof. Let fM equal / on the support [—M,M] k and zero outside. For the choice of contours 7 in the 
definition of there exists C \, C 2 such that C\ < < C 2 holds for all z £ 7. This along with 

the exponential growth bounds on / implies that {fF q ' v /m)(<?) is uniformly convergent as the Zi vary 
along 7. This, along with the fact that J’ q,u J :q,u fM = fM implies the desired result. □ 

Definition A.5. Let W]J iax be the space of all functions / : WY fc —>■ C such that J q ' v J :q ^f = f (and 
likewise for Wf ;max ). 

The following corollary is a consequence of Proposition IA.41 
Corollary A.6 . For c > 0 small enough and C > 0 large enough, H^ xp( - C ^ C (and likewise for 

w£ ). 

1 ;max/ 


Appendix B. Vertex weights for J = 1,2,3 


B.l. J = 1 vertex weights. 


When g 



B.2. J = 2 vertex weights. 
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Again, note the automatic vanishing of suitable probabilities triggered by factors 1 — q 9 and q — q g . 


B.3. J = 3 vertex weights. For J = 3, there are 16 vertex types, so we will no longer draw the arrow 
configurations. We will also omit the common denominator (1 + a)(l + qa)( 1 + q 2 a) which is present 
in all probabilities. The table of J = 3 vertex weights is the following: 



h = 0 

32 = 1 

32=2 

CO 

II 

CM 

h = 0 

(1 + aq 9 )(l + aq 9+1 ) 
x (1 + aq 9+2 ) 

a{l + q + q 2 )(l + aq 9 ) 
x (1 + aq 9+l ){ 1 — q 9 ) 

a 2 {l + q + q 2 ) 

X (1 + aq 9 ) 

X (1 — q 9 )(q — q 9 ) 

a 3 ( 1 — (/ 9 ) 

X (q — q 9 )(q 2 — q 9 ) 

31 = 1 

(1 — vq 9 )(l + aq 9+1 ) 
x (1 + aq 9+2 ) 

(1 + aq 9+l ) 

X ^a(l + q + q 2 ) 

— aq 9+1 ( 1 + q — aq) 
+ z^g 9 (l — a — aq) 

+ avq 29 (l + q + q 2 ) S j 

a 2 q (l + q + q 2 ) 

+ aq 9 (aq — q — 1) 

X (—z/ + cm/ 2 + aq) 

+ aq 2g (a + aq — 1) 

X [u — aq 2 + vq) 

+ a 2 vq ig (l + q + q 2 ) 

a 2 (l — q 9 )(q — q 9 ) 

X (aq 2 + vq 9 ) 

3i = 2 

(1 - vq 9 )(l - vq 9+1 ) 
x (1 + aq 9+2 ) 

(1 — vq 9 ) 

X ^a(l + q + q 2 ) 

+ aq 9+2 (a + qa— 1) 
+ vq 9 ( 1 + q — qa) 

+ cwg 2si+1 

X (1 + g + q 2 )^j 

(aq + vq 9 ) 

X ^a(l + q + q 2 ) 

— aq 9+1 (1 + q — qa) 
+ vq 9 ( 1 — a — qa) 

+ avq 2a (l + q + q 2 ) S j 

a( 1 — q 9 )(aq + vq 9 ) 

X (aq 2 + vq 9 ) 

3i = 3 

(1 - vq 9 )(l - vq 9+L ) 

X (1 - uq 9+2 ) 

(1 + q + q 2 )(a + vq 9 ) 
x (1 — z^g 9 )( 1 — ^( 7 9+1 ) 

(1 + q + q 2 )(a + vq 9 ) 

X (aq + vq 9 )( 1 — vq 9 ) 

(a + vq 9 )(aq + vq 9 ) 

X (aq 2 + vq 9 ) 
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Appendix C. Yang-Baxter equation 


Denote, for this section, R^’al{k\.k 2 '-k\. k' 2 ) = R^}(rn,k\\£. k\ ) R,ijJ(£, k‘ 2 ] n, ki 2 ) where £ 


k[ = n + k' 2 - k 2 . Let R^’al(k\,k 2 ; k[,k! 2 ) = R^’ai(k[, 




k '2 i k\, 


^2) and define the matrix 


m + k\ — 


/1 

0 

0 

q(u\-u 2 ) 

su\u 2 (u\—qu 2 ) 

n 

1-q 


su 2 {u\—qu 2 ) 

1 0 

0 


0 

1 -Q 

sui(ui—qu 2 ) 

u\—u 2 

su\u 2 (u\—qu 2 ) 

0 


0 \ 
0 

0 

1 

S 2 u'fu'2 ) 


under the association on. = — sui and v = s 2 . Then the Yang-Baxter equation amounts to the fact that 
Roi^’al and Rq'j 'oI (note the interchange of indices in or, 0:2) are similar with respect to Y: 


r>(m,n)y _ y T>(m,n) 

- tly Oi2,Q!i ±k 'Oli,OL2’ 

This can be derived, in light of Remark 12.21 from |Borl4l Proposition 2.5], which itself is just a restate¬ 
ment of the standard six-vertex Yang-Baxter equation. 
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